Mon. Not. R. Astron. Soc. 0Q0.mi2T](2Q10) Printed 20 October 201 1 (MN KT^ style file v2.2) 



The velocity modulation of galaxy properties in and near clusters: 
quantifying the decrease in star formation in backsplash galaxies 

^ ; Smriti Mahajan^*, Gary A. Mamon^'^ and Somak Raychaudhury^ 

»' ' 1: School of Physics and Astronomy, Universitv of Birmingham, Birmingham B15 ITT, United Kingdom 
2: Institut d'Astrophysique de Paris (lAP), 98-bis Blvd Arago, F-75014 Paris, France 

3: Astrophysics & BIPAC, Department of Physics, University of Oxford, Keble Rd, Oxford OXI 3RH, United Kingdom 



ABSTRACT 

Ongoing/recent star formation in galaxies is known to increase with increasing projected dis- 
tance from the centre of a cluster out to several times its virial radius (i?v)- Using a complete 
sample {Mr ^ —20.5, 0.02 ^ z ^ 0.15) of galaxies in and around 268 clusters from the 
SDSS DR4, we investigate how, at a given projected radius from the cluster centre, the stellar 
mass and star formation properties of a galaxy depend on its absolute line-of-sight velocity 
in the cluster rest frame, jtiLosI- We find that for projected radii R < 0.5 Rv, the fraction of 
high mass non-BCG galaxies increases towards the centre for low |wlos|; which may be the 
consequence of the faster orbital decay of massive galaxies by dynamical friction. At a given 
projected radius, the fraction of Galaxies with Ongoing or Recent (<l-3 Gyr) Efficient Star 
Formation (GORES, with EW(H(5) > 2 A & D„4000 > 1.5) is sUghtly but significantly lower 
for low |f LOS I galaxies than for their high velocity counterparts. We study these observational 
trends with the help of a dark matter (DM) cosmological simulation. We classify DM parti- 
cles as virial, infall, and backsplash according to their present positions in (r, Vr) radial phase 
space and measure the frequencies of each class in cells of {R, |f los I) projected phase space. 
As expected, the virial class dominates at projected radii R < R^, while the infall particles 
dominate outside, especially at high |wlos|- However, the backsplash particles account for at 
least one-third (half) of all particles at projected radii slightly greater than the virial radius 
and li'LOsI < <^v (I'I'losI ^cr„). The deprojection of the GORES fraction leads to a saturated 
linear increase with radius. We fit simple models of the fraction of GORES as a function of 
class only or class and distance to the cluster centre (as in our deprojected fraction). While 
GORES account for 18 ± 1% of all galaxies within the virial cylinder, in our best-fit model, 
they account for 13 ± 1% of galaxies within the virial sphere, 11 ± 1% of the virial popula- 
tion, 34 ± 1% of the distant (for projected radii R < 2 R^) infall population and 19 ± 4% of 
the backsplash galaxies. Also, 44 ± 2% of the GORES within the virial cylinder are outside 
the virial sphere. These fractions are very robust to the precise good-fitting model and to our 
scheme for assigning simulation particle classes according to their positions in radial phase 
space (except for two of our models, where the fraction of GORES reaches 27 ± 4%). Given 
the 1-3 Gyr lookback time of our GORES indicators, these results suggest that star formation 
in a galaxy is almost completely quenched in a single passage through the cluster. 

Key words: Galaxies: clusters: general; galaxies: evolution; galaxies: starburst; galaxies: star 
formation; galaxies: kinematics and dynamics 



1 INTRODUCTION 

The wide variety in star formation properties of galaxies is 
an inevitable consequence of the varied environments they 
inhabit. Physical propertie s of galaxies such as morphology 



Balogh et alj|2004bh . luminosity (e.g. 



I I 19981) . individual spectral indices used 
as star formation rate (SFR) tracers, such as the equivalent width 
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(EW) of [Oil] or Ha (e.g. lSalogh et alj2004al : lHaines et al.ll2006h . 

are known to vary systematically with environment. As they leave 
their sparsely populated base in the field for the high density re- 
gions of galaxy clusters, galaxies see their dominant stellar popu- 
lation become older, emission lines disappear from their spectra, 
complex spirals are transformed into more uniform spheroids, and 
the SFR decreases by 1-2 orders of magnitude. 

It is thus generally believed that the cluster environment is 
hostile to star formation. For example, the denser the environ- 
ment, the stronger will be the stripping of interstellar gas either 
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by the ram pressure of the intra-c luster ga s llGunn&Gotllll972h 
or by the tidal field of the cluster iMerritll 1983 ), leading to star- 
vation for subsequent star formation ( iLarson. Tinslev. & Caldwelll 
Il98d) . So a galaxy that entered the cluster virial radius ~ 1 
Gyr ago, and has since crossed the core of the cluster, should 
have lost its gas via either of these two effects during it s pas- 
sage through the cluster core (see for instance IVollmej |2009| 
for the impact of ram-pressure stripping in the Virgo cluster). 
Moreover, in the intermediate density environments on the clus- 
ter outskirts , mech anisms like galaxy-galaxy interactions (e.g. 



evolution (e.^. [ Porter & Ra v chaud hur^ 



falling galaxies are expected to have distinctive properties com- 
pared to the passively evolving virialised population inside the 
cluster. 

However, one can define an intermediate population of galax- 
ies that have travelled through the cluster core (perhaps more than 
once) and have not yet had time to mix with the virialised popu- 
lation. These backsplasl^ galaxies would have felt the first effects 
of the cluster environment, but may not have evolved all the way 
to attain the passive proper ties of the virialised populat i on. In their 
simulations of six clusters, iBalogh. Navarro. & MorrisI ( |2000|) find 
that 54% ± 20% of dark matter (DM) particles at 1 - 2 R200 from 
the cluster c entre have been i nside the virial radius (i?20o) at some 
earlier time, et alj ( l2005h estimate similar fractions, consider- 
ing a virial radius equivalent to -Rioo, and further add that 90% of 
such p articles go as close to the cluster centre as 50% of the virial 
radius. iMamon et aL used both analytical methods and cos- 

mological simulations to conclude that DM particles travel through 
the core of the cluster out to 1-2.5 J?ioo on the other side. 

The backsplash galaxies may contribute to the observation 
that the radial trends of st ar formation diagno stics extend to at 
least twice the virial radius tealogh et ai]|2000l) from cluster cen- 
tres. Hence, the intermediate environments found in inter-cluster 
filaments and on the outskirts of galaxy cluste rs play a key role 



in defining the e volutionary history of galaxies llRines et"al] 20031 : 



Gill et al.ll2003: IPorter & Ravchaudhurvl I200I I Porter et al.l 
Mahaian et al.l201 11201 ih 



20081 : 



Most results indicate that the properties of galaxies evolve as 
they encounter denser environments, i.e. progress from the field to- 
wards the cluster. This could be due to the prevalence of smaller 
clusters (galaxy groups) as one approaches rich clusters, as the 
tidal field and possibly some mild form of ram pressure stripping 
might play a similar role, as they do in more massive clusters, to 
deplete the gas reservoirs of galaxies. Alternatively, it is possible 
that the gas reserves of galaxies may be protected by the pressure 
of the surrounding intra-group medium (IGM) from being tidally 
perturbed, preventing the star formation in the galaxy from being 
quenched, even after spending > 1 Gyr in crossing the cluster. In 
such a case, an increase in tides generated during the group-cluster 
merger may instead re sult in an enhancement of the SFR of group 
galaxies ( lBekkilll999h . Although such a group pre-processing sce- 
nario is yet to be confirmed, som e evidence in i ts favour is found in 
recent observational studies (e.g. lOemler et al.ll2009l) . 

The main difference between the infalling and backsplash 
galaxies is that the former have larger apocentres, hence larger ra- 
dial (3D) velocities. Hence, one expects that the backsplash galax- 



ies will have lower absolute line-of-sight (hereafter, LOS) veloci- 
tie|f|than the infalling galaxies. One is therefore motivated to study 
the modulation of the radial trends of galaxy properties, with abso- 
lute LOS velocity. 

In the central regions of clusters, the relative velocities of 
gala xies have been found to be correlated with their luminosity 
(e.g.lChincarini & Roodlll977l : IStrublell97^:lBothun & Schomben 



1990"; 'Biviano et al.' '1992'). morphology de V aucouleurd I Il96ll 
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galaxies are found to have lower velocity dispersion relative to 
their late-type, star-forming counterparts. It has been shown that 
this velocity segregation of mass (or luminosity) in cluster galaxies 
is neithe r induced by morpho logical segregation nor limited to cD 
galaxies iBiviano et al.il 9921 ). 

There have been some efforts to study the velocity segrega- 
tion in galaxy properties on the outskirts of galaxy clus ters (e.g. 
iMohr et alJI 19961: iRines et alj|2005l : [Pimbblet et al]|2006l) . In par- 
ticular, Rines et al. 1 20051) found no significant difference in the 



velocity distribution of the emission-line and absorption-line galax- 
ies with projected radii between 1 and 2 -R200 in a sample of eight 
nearby x-ray bright clusters. However, emission-line galaxies out- 
side the virial radius appear to have a significantly lar ger velocity 
dispersion than their absorption-line counterparts (but lRines et al.l 
do not provide any quantitative e stimate of the significance). More- 
over, from fig. 6 of iRines et all the mean velocity dispersion pro- 
file of emission-line galaxies appears significantly steeper than that 
of non -emission-li ne galaxies (this trend was not noted by the au- 
thors). |Rines|eral] concluded that star-forming galaxies in the infall 
regions are not field interlopers, which lead them to reject a simple 
model where the truncation of star formation at low projected ra- 
dius was solely caused by the efficient quenching of star formation 
in (backsplash) galaxies that have crossed through the cluster. 

In this paper, using the very large Sloan Digital Sky Survey 
(SDSS), we have collected sufficient data around hundreds of clus- 
ters to analyse the differences in the distribution in projected phase 
space (projected radius R and absolute LOS velocity |wlos|) of 
star-forming and passive galaxies at much finer resolution than pre- 
viously performed. We look for and find a significant velocity mod- 
ulation of the variations in stellar mass and star-formation proper- 
ties with projected radius. We also build simple dynamical models 
involving the virial, infalling and backsplash populations of galax- 
ies to explain our results, calibrating our models and the projection 
effects using a cosmological simulation. 

The organisation of this paper is as follows: in the following 
section we describe the observational dataset and briefly describe 
our star formation diagnostics. In ij3]we present the analysis of the 
SDSS galaxies around clusters, focusing on the velocity modula- 
tion of radial trends in stellar mass and star formation efficiency. 
In ij4]we model our results using a cosmological simulation de- 
scribed in t]4. 1 I to classify particles into virial, infall and backsplash 
populations ( t]4.2| l. In i|4.3l we fit simple models of the variation 
with physical radius of the fraction of Galaxies with Ongoing or 
Recent Efficient Star formation (hereafter, GORES) for each of the 
three populations to the observed fraction of GORES as a function 
of projected radius and LOS velocity. We discuss and summarise 
our results and their implications in ^ Throughout this paper we 



The term backsplash galaxies was first coined bvlGiiT Knebe. & GibsonI 
i2005l) . 



^ Our absolute LOS velocities are in the cluster frame. 
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assume a concordance cosmology with Sin 
h = Q.7. 



■- 0.3, Q.A = 0.7, and 



2 DATA 

2.1 The cluster catalogue 

We sele cted our ga l axy cl usters from the 'group' catalogue com- 
piled bv lYang et allbOOTi) from the Fourth Data Release (DR4) of 
the SDSS. Their iterative algorithm identifies first order groups us- 
ing a Friends-of-Friends (FoF) algorithm in redshift space, with a 
small linking length in projected space and a large one along the di- 
rection of the LOS. This algorithm thus selects potential groups at a 
high surface density threshold, so that groups embedded within the 
virial radii of more massive ones will be extracted as two (or more) 
separate entities. All galaxies left out by this grouping method are 
treated as additional potential groups (of single galaxies). 

The raw luminosities and (luminosity-weighted) centres of 
the clusters are measured and then the luminosities (corrected 
for incompleteness), masses (assuming a given mass-to-light ra- 
tio) and velocity dispersions (from a mass / velocity dispersion 
relation) are derived. The radii -Riso, i.e. of spheres that are 180 
times denser than the critical density of the Universe are de- 
rived from the masses. The distribution of galaxies around each 
group centre in projected phase space are examined to determine 
group membership, based upo n the assumption of a projected 
iNavarro. Frenk. & White! ( Il996l . hereafter NFW) surface density 
profile and a Gaussian distribution of LOS velocities with disper- 
sion independent of projected radius. The algorithm is iterated with 
these new group assignments, where the mass-to-light ratio used for 
estimating the new mass is a function of luminosity derived from 
the previo us iteration. The ite ration continues until the membership 
converges. lYang et al.l J2007h have tested their algorithm on mock 
galaxy catalogues constructed in real space and converted to red- 
shift space. 

In the present work, we define the galaxy clusters with the 
following criteria: 

(i) redshift range: 0.02 ^ 2 ^ 0.12; 

(ii) membership: at least 15 galaxies; 

(iii) minimum halo mass (within i?i8o): IO^^Mq. 

All of our chosen clusters have at least 12 galaxies (Mr ^ —20.5) 
with SDSS spectra within i?i8o of each cluster centre. For our cho- 
sen cosmology, the virial radius, within which the cluster may be 
thought to be in dynamical equilibrium, is commonly believed to 
be very close to -Rioo where the density is 100 times critical. In this 
paper, we adopt -Rioo for the virial radius, and more precisely 



Rv — 1.3 Ri 



(1) 



which, for NFW models with concentrations between 3 and 10, 
matches iiioo to an accuracy of 2%. This gives us a sample of 
268 clusters in the SDSS DR4 footprint with median Miso ~ 
1O^*'^M0, corresponding to a median virial radius -Rioo = 
1.5 Mpc (using Eq.lTJ. 

The range of 1.5 dex in cluster masses should not be a con- 
cern given that the properties of non-central (satellite) galaxies 
depend little on the mass of the group or cluster they live in 
jvan den Bosch et al.ll2008l) . 



2.2 Galaxy data 

The observed photometric and spectroscopic data for 

[ alaxies in this work are taken from the SDSS DR4 
Adelman-McCarthv & the SDSS coUaboratiorJ |2006|) . We 
use the galaxy magnitudes, fc-corrections and the corresponding 
galactic extinction values in the SDSS g and r-bands from the New 
York University Va lue Added Galaxy Catalogue (NYU-VAGC; 
lBlantonetaLll2005l) . We fc-correct these magnitudes to 2: = 0.1, 
which is the median redshift of the sample. We select galaxies that 
meet the following criteria: 

(i) redshift range 0.02 2 0.15; 

(ii) absolute magnitude: Air ^ —20.5; 

(iii) half-light angular radius: Oeft ^ 5"; 

(iv) projected separation with cluster: R < 2R^\ 

(v) velocity relative to cluster: |wLOS-"Los°'^i ^ 3 0-^'"'''°''. 

Our luminosity limit roughly c orresponds to L ^ 0.5 I/* (given 
that M^:l = -20.44 4-5 log fc. lBlanton et"al]|2003h . which corre- 
sponds to the apparent magnitude limit (r ^ 17.77) for the SDSS 
spectroscopic catalogue at 2 = 0.1. Our galaxy redshifts are chosen 
in an interval that is slightly larger than for the clusters, so as to 
avoid missing galaxies lying behind (or in front) of the most (least) 
distant clusters. The angular size criterion ensures that the spectrum 
obtained from the 3" SDSS fibre is a fair representation of the light 
from the galaxy, and does not just correspond to the passive central 
bulge of the galaxy. 

This gives us a sample of 19,904 galaxies in or near 268 clus- 
ters covered out to 2 virial radii in projection and 3 times the cluster 
velocity dispersions along the LOS. 

The cluster galaxies are compared with a sample of 21,100 
field galaxies, randomly chosen to have the same redshift and lu- 
minosity distributions as those of the cluster galaxies. We require 
the field galaxies to lie at least 10 Mpc away from the c entres of all 
group s with N ^ 5 and halo mass ^ 10^^'^ M© in the I Yang et al.l 
( I2OO7I) catalogue (this includes clusters). Since the galaxy prop- 
erties degend_on_their_distance from the nearest large structure 
(e.g. lMahaianet alJl20Ill) out to ~ 3i?2oo (< 6 /i^ Mpc; e.g. 
iRines et alj2003l ; lMahaian et al.l201ll) . a linear distance of 10 Mpc 
(corresponding to typically 6.7 R^) from all groups and clusters en- 
sures that our field sample is not contaminated by the effect of any 
groups or clusters. 



2.3 Star formation diagnostics 

In this paper, as our principal diagnostic, we use the specific star 
formation rat e (SFR/M*) of a galaxy, a s estimated from the optical 
spectrum, by iBrinchmann et alj l l2004h for all galaxies belonging 
to the spectroscopic catalogue of SDSS DR4. From this and a few 
other spectral indices, namely the emission-corrected USa (here- 
after, US) EW and the 4000 A break, estimated in a related work 
jKauffmann et al.ll2003l) . we compare the trends seen in the cur- 
rent star formation activity and star formation histories (SFHs) of 
diffe rent types of galax ies. 

IBrinchmann et al.l have construct ed a model grid of ~ 
2 X lO'^ spectral synthe sis models ( iBruzual & Chariot! !l993! ; 
!Charlot & Longhetti!l200ll) over four main parameters spanning a 
range of -1< log Z/Zq < 0.6 in metallicity (Z), -4.0 < log f/ < 
—2.0 in ionisation parameter (U), 0.01 < tv < 4.0 in total dust 
extinction (tv) and 0.1 < ^ < 0.5 in dust-to-metal ratio (^). 
Dust attenuation is al so treated with a very sophisticated model of 
!Charlot & Falil l !2000l) . which provides a consistent model for ul- 
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0.5<r/R^gl.O 1.0<r/R^gl.5 1.5<r/R^S2.0 




•g 9 10 11 9 10 11 9 10 11 9 10 11 
log m'/Mq log m'/Mq log m'/Mq log m'/Mq 

Figure 1. The distribution of stellar masses for our galaxy samples, in four bins of projected radius in the four panels from left to right. The solid black, dashed 
red and dotted blue distributions represent absolute LOS velocity bins Ii'losI/o'u = 0-1, 1-2 and 2-3 respectively. The thin grey distribution denotes the 
field. The bottom panels show the difference between the binned cluster distributions and that of the field. The error bars are Poisson. In the flux limited SDSS, 
galaxies in the local Universe have their M* distribution centred around A/ * = 3 X lO^" M0 teauffmann 6710120031) . A significant fraction of galaxies 
with A/* ^ 10^^ Mq are found only in the core of galaxy clusters. Interestingly, the low absolute LOS velocity galaxies in the cores are the ones which are 
more massive, while the high velocity galaxies mostly inhabit the low mass end. 



traviolet (UV) to infrared (IR) emission. Tiiey then go on to derive 
tlie values for the SFR and SFR/M* using a Bayesian technique, 
considering the entire spectrum of each galaxy. The derived pa- 
rameters (SFR and SFR/Af*) are then corrected for aperture biases 
using photometric colours. Their analysis of line ratios is insensi- 
tive to stellar ages, SFH and the relative attenuation by dust in the 
gas clouds and the interstellar medium. Out of the three statistical 
estimates (mean, median & mode) for the probability distribution 
function of SFR/ A/* derived for each galaxy, we used the median 
of the distribution in this paper, since it is independent of binning. 

The spectrum of a galaxy that has experienced a burst of star 
formation in the recent past (1-3 Gyr ago) is expected to show 
a strong absorption of US in the optical spectrum. This is be- 
cause A type stars, which significantly contribute to the HS ab- 
sorption, have a lifetime of < 1.5 Gyr, but a burst of star forma- 
tion can produce enough ionising radiation to fill in the US ab- 
sorption. However, as the episode of star formation ends, the US 
absorption line becomes more transparent. The EW(H(5 ) has been 
extensively used in the literature to measure the mean stellar age 
of galaxies, and the elapsed time since the last m ajor burst (e.g. 
IWorthev & Ottaviani|[l997l : iKauffmann et alj|2003l) . At the other 
extreme is the 4000 A break (D„4000), which is the strongest dis- 
continuity occurring in the optical spectrum of a galaxy, from the 
accumulation of absorption lines of mainly ionised metals. As the 
opacity increases with decreasing stellar temperature, the value of 
D„4000 is largest for old and metal-rich stellar populations. Stel- 
lar population synthesis models indicate that D„4000 is mostly 
sensitive to age, except for sub-solar metallicity populations older 
than 1 Gyr, where it is also strongly dependent on metallicity 
jKauffmann et alj2003h . 

Interestingly, the variation of H<5 with D„4000 is not sensi- 



tive to metallicity or dust extinction, thus providing a good mea- 
sure of SFH, such that for a given D„4000, stronger H5 im- 
plies th at the galaxy has exper ienced a recent burst in star for- 
mation i Kauffmann et al.|[2003h . In particular, the right panel of 
IKauffmann et all 's fig. 2 indicates that the combination of US > 
2 A and D„4000 < 1.5 leads to stellar populations younger than 
1 to 3 Gyr for solar and 0.2 solar metallicity populations, respec- 
tively. 



3 VELOCITY MODULATION OF GALAXY 
PROPERTIES 

The key issue that we are trying to address in this paper is whether 
it is possible to distinguish the galaxies that have passed through 
the cluster core from those that are falling into it for the first time, 
based on observable parameters. DM simulations show that the 
backsplash population is most likely to be seen between 1 — 1.5i?v 
and with velocities l ower than the m ean velocity of the cluster (e.g. 
iMamon et ai]|2004 loill et alj|2005l) . So, we stack together all the 
galaxies found in clusters within 2i?v and |Aulos1 < 3(t„, into 
four bins of projected radius of width 0.5-Rv, and into three bins of 
absolute LOS velocity: | Aulos |/o"i> <1, 1-2 and 2-3 respectively. 
Below we describe the features of these distributions for various 
galaxy properties with respect to one another and the field distribu- 
tion. 



3.1 Stellar mass 

In Fig.[T]we plot the distribution of stellar masses for the 12 classes 
(3 scaled velocity bins in each of the 4 scaled radial bins) of galax- 
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ies as described above. Due to the flux limit of SDSS, the stellar 
mass distribution of galaxies in Fig.[T]appears to be centred around 
M* = 3x10^" Mq in the local Universe (also see lKauffmann et aP 
l2003h . The most massive galaxies (log Af */Mq ^11) are mostly 
confined to the cores of the clusters. However, even in the core 
there seems to be a correlation between the relative LOS veloc- 
ity of a galaxy with respect to the cluster and its stellar mass. The 
Kolmogorov-Smimo\[f| (K-S henceforth) statistics show that the 
most massive galaxies in clusters have low velocities with respect 
to the cluster, and the high velocity galaxies show an affinity for 
the field distribution, which is skewed towards the low mass end. 
The detailed statistical comparison between different distributions 
is shown in Table lAll in Appendix. 

When the brightest cluster galaxies (BCGs) are excluded, 
there is no mass segregat ion in clusters relative to the field 
jvon der Linden et alj 2010l) . To explore the inter-dependence of 
M*, Rv and LOS velocity of galaxies in clusters, we repeated 
o ur analysis by excluding the BCGs from the clusters. Contrary 
to Ivon der Linden et alj | |2010|) . we find that the different distribu- 
tions still show similar trends: in the inner regions of clusters the 
most massive galaxies have lowest absolute LOS velocities (see Ta- 
blelA2t. 

To analyse mass segregation in a different way, in Fig.|2]we 
show the fraction of high-stellar mass galaxies as a function of 
projected radius in three bins of absolute LOS velocity, discard- 
ing BCGs. Fig.[2]confirms our impression that mass segregation is 
present, but only at R < 0.5 -Rv, and only for the low- and inter- 
mediate |«LOs| galaxies. 

The difference in low and high mass velocity distributions in 
the projected cluster core is not simply caused by the BCGs, which 
are ex cluded from Fig.[2](see also Table lA2T ). While Isiviano et al] 
( Il992h also found that high mass non-BCG galaxies had lower ab- 
solute LOS velocities, we find that this is limited to galaxies within 
0.5Rv In other words, we find that the high absolute LOS veloc- 
ity galaxies show no mass segregation, in cont rast with the global 



luminosity segregation fou nd by others (e.g. lAdami et alj 



Ivon der Linden etai]|2010l) . when no cuts are made on absolute 
LOS velocity (which we also find, since the high absolute LOS 
velocity galaxies are much rarer than their intermediate and low 
velocity counterparts). We discuss mass segregation amongst clus- 
ter galaxies further in i]5.1l 



3.2 Galaxy colours 

Broadband photometric colours are widely used as proxies for the 
star formation activity of galaxies. In this paper, we use the 'colour 
offset', which is the difference between the — J")"'^ (extinction- 
corrected) colour of a galaxy and that of the red sequence for a 
galaxy of the same magnitude, where the red sequence is fitted by 
stacking all the cluster galaxies together. In Fig. [3] we plot the dis- 
tribution of the colour offset for all 12 classes of galaxies. 

The low velocity galaxies in cluster cores are found to have 
the reddest colours, indicative of their passive evolution. But inter- 
estingly, even within the cores, the difference between the colours 
of high and low velocity galaxies indicates statistically significant 
different distributions (see Table lASV Out to 2Rv, the distribution 



^ The null hypothesis for the K-S test is that both the distributions be- 
ing compared ai'e drawn from the same parent distribution, i.e. if the two 
distributions are identical, the probability (p) of the null hypothesis being 
satisfied is unity. 




Figure 2. The fraction of high stellar mass galaxies (M* 3 X 10^° Mq, 
with brightest cluster galaxies removed), plotted as a function of scaled ra- 
dius from the cluster centre. The low (0 — 1 (T„) and high (> 2 (t„) absolute 
velocities (shown as black circles and blue diamonds respectively), and the 
±1 cr range (from binomial statistics) are shown as black dash-dotted lines, 
red dashed lines, and blue error bars, for the low, intermediate (1-2 a^ ) and 
high LOS absolute velocities, respectively. The magenta arrow indicates the 
field value. 



of the colour offset for the field galaxies is very different from that 
of the galaxies residing in and around clusters. This might sug- 
gest that the broadband colours could be very easily modulated 
by slight changes in the local and/or global environment, and are 
hence not very well suited f or the study of evolutionary trends 



amongst galaxi es (see also Mahaian. Haine s. & Ravchaudhurvl 



[2010; Ma haian. Ravchaudhurv. & Pimbblet .201 ih . As discussed 
below, on the outskirts of clusters (l-1.5i?v) the trends seen for 
broadband colours are consistent with those seen for spectroscopic 
parameters. 



3.3 Diagnostics of specific star formation rate and star 
formation iiistory 

It has been shown that the SFR of galaxies is a function of their 



stellar mass (e.g. lNoeske et alj2007bl) . So in this work we prefer to 
use the specific star formation rate (SFR/M*) to quantify the cur- 
rent star formation activity in a gala xy. In Fig. [4] we plot the d istri- 
butions of SFR/ A/*, as estimated bv lBrinchmann et al] ( l2004h . for 
all the galaxies in the spectroscopic galaxy catalogue of the SDSS 
DR4. 

As expected, the cluster galaxies within 0.5 7?v have their log 
SFR/A/* distributions heavily skewed towards the lower end, in- 
dicating passive evolution. Interestingly though, the high end of 
the values of SFR/Af * in the cluster core, and almost at all radii 
out to 1.5 J?v, mostly belongs to galaxies with high velocities. K- 
S statistics show the difference between the lowest and the high- 
est velocity galaxy distributions to be statistically significant in the 
core (< 0.5 7?v) and on the cluster outskirts (1 — 1.5 R^; see Ta- 
ble IA4t . The field distribution is significantly different from all 
other classes of galaxies, becoming indistinguishable from the high 
velocity galaxies only in the outermost radial bin. 
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Figure 3. Same as Fig.[T]but for the "colour offset" of cluster galaxies, i.e., the difference between the (g — r)"'^ colour of a galaxy and the Red Sequence 
colour for a galaxy of the same magnitude. Galaxies bluer than the fitted red sequence (see text) have positive values of the offset on the x-axis. The build-up 
of the blue tail of the distribution outwards from the cluster core is clearly seen. Interestingly, even in the outermost radial bin (1.5-2JJv), there are more 
red sequence galaxies than the field. This shows that physical processes, dependent on environment, efficiently modify properties of galaxies in regions with 
densities in between the two extremes characterising cluster cores and the field. 
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Figure 4. Same as Fig.[T]but for SFR/M*. As expected, the cluster distributions are skewed towards the low log SFR/Af* end in the core region, but only a 
small tail of galaxies with low log SFR/M* is seen in the outer radial bins. 
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Figure 5. Same as Fig.[T]but for EW(H5). This plot complements the results in Fig.|5] The field distribution shows a clear bimodality around 2 A. The galaxies 
in the core of the cluster, most of which are evolving passively, do not show absoiption in H5. It is interesting to note that the high velocity galaxies in the 
innermost bins, which are likely to be the newest entrants in the core, are also the ones which are most likely to have experienced a burst of star formation as 
indicated by relatively high EW(H<5). This trend in velocity classes continues out to 2Rv, where the difference between the field and the low velocity galaxies 
is minimum. The difference between the highest (blue) and the lowest velocity (black) galaxies is statistically significant in all radial bins out to 1.5 Rv 
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Figure 6. Same as Fig. \T\ but for Dn4000. The field distribution is bimodal with a minimum at D„4000 
~ 1.5. The passively evolving galaxies in the cluster core are old and metal-rich, hence their distribution peaks at the high end of the D„4000. How- 
ever, the high velocity galaxies in the cluster core (dotted blue line) have a distribution similar to that of field galaxies. The trend continues out to 1.5 Rv, 
beyond which the distributions almost merge with the field. The K-S statistics show the distributions of the high and low velocity galaxies to be statistically 
different at all radial distances within l.Siiv. 
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r/R^<=0.5 
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Figure 7. The distribution of galaxies having the highest (\vi^os\/ <Jv = 2 — 3; blue points) and the lowest (Ii^losI < ""d; black contours) absolute LOS 
velocity relative to the cluster, in bins of increasing (scaled) radius (clockwise from top left). The innermost contour represents a density of 100, decreasing by 
10 at every subsequent level. The fractions in the top right corner of each panel show the ratio of the numbers in the two classes of galaxies shown, and the 
scatter in it. Barring the innermost radial bin, the distribution of high velocity galaxies remains the same, although the distribution of the virialised galaxies 
tends to become bimodal with increasing (scaled) radius. 



Figs.|5]and|6]show the distributions of EW(H5f| and D„4000 
for the galaxies in 12 different cells of projected phase space 
(three scaled absolute LOS velocities times four scaled projected 
radii), and compare them with the distribution of field galaxies. The 
EW(H5) distribution of the cluster galaxies is strikingly different 
from that of the field galaxies, the latter showing a bimodality at 
EW(H5) ~ 2 A (also see Table lASl l. In the cluster cores, the galax- 
ies with highest velocities, which are most likely to have fallen in 
recently (Table 1), are the ones that also show H(5 in absorption 
(EW(H<5) > 0). This suggests that their likelihood of experiencing 
a burst of star formation within the last ~ 1 — 3 Gyr is relatively 
higher than their low velocity counterparts. The K-S test shows that 
the difference between the EW(H5) distributions of the highest and 
the lowest absolute LOS velocity galaxies is statistically significant 
at all radii within L5 i?v (TablelASll. 

The amplitudes of the 4000 A break in SDSS galaxy spectra 
show very different distributions for the highest and lowest classes 
of |wlos|, both in the core (< 0.5 i?v) and on the periphery (1 — 
1.5 7?v) of the clusters. However, the distributions of D„4000 in 



* Note that the negative values of EW(H(5) do not imply H5 line emission 
in the passive galaxies, because the EW(H(5) values are measured from the 
pseudo-continuum bands, and so do not always represent the 'true' flux 
in the H<5 line. The pseudo-continuum bands are required because in the 
spectra of passive galaxies, the position of several metal lines may coin- 
cide with the H(5 line, making it virtually i mpossible to measure the 't rue' 
flux (G. Kauffmann, p riv. comm.; see also lWorthev & Ottavianiiri997l and 
iKauffmann et alj2003h . 



the 0.5 — 1 i?v bin are only very marginally dissimilar, with ~ 10% 
probability of the null hypothesis being satisfied (see Table |A6t . 

We now take advantage of the fact that together and 
Dn40 00 yield a good constraint on the SFHs ( Kauf fmann et al.l 
l2003h of galaxies. In Fig.|7]we show the SFHs of galaxies with dif- 
ferent absolute LOS velocities in bins of increasing cluster-centric 
radius. For clarity, we only plot the two extreme absolute LOS ve- 
locity classes (|t;LOs| < and |ulos|/'7"u = 2 — 3). The distribu- 
tion of the low \v-Los I galaxies, in this plot of EW(H(5) vs D„4000, 
is dominated by the passive sequence (lower right) at low projected 
radii, and becomes increasingly bimodal with increasing cluster- 
centric radius. The same evolution occurs for the high |i'los | galax- 
ies, but is less bimodal at large projected radius. On the other hand, 
the high absolute LOS velocity galaxies span the entire available 
parameter space e qually at all distan ces. According to the stellar 
population study of lKauffmann et akl (see their fig. 3), this is an in- 
dication that the SFH is more continuous in the low |ulos | galaxies 
and star formation occurs more in bursts in the high velocity galax- 
ies. We discuss these trends in the context of the relative fraction 
of the infalling and backsplash galaxies, and the fraction of star- 
forming galaxies as a function of cluster-centric distance in i|5] 

We now refer to galaxies with EW(H5) > 2 and D„4000 < 
1.5 as Galaxies with Ongoing or Recent Efficient Star formation 
(GORES), following the assumption that a galaxy that experienced 
efficient star formation < 1 — 3 Gyr ago is likely to show these 
features prominently in its spectrum. Our choice of EW(H5) and 
Dn4000 thresholds are based on the bimodality seen in the field 
distributions (Figs.[5]and[6j. 

Fig. [8] shows how the fraction of GORES, /gores, varies 
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Figure 8. The fraction of Galaxies witli Ongoing or Recent Efficient Star 
Formation (GORES, EW(H(5) > 2 A and D„4000 < 1.5) versus projected 
radius in bins of absolute LOS velocity: — 1 o-„ {black triangles), 1 — 2 Cj, 
{red squares), and > 2 (t„ {blue pentagons), slightly shifted along the x- 
axis for clarity. The magenta arrow indicates the GORES fraction in the 
field {R > 3i?v). 



for subs equent epochs (i.e., at 7.5 h^^ kpc). DM halos were iden- 
tified bv lBorgani et alj at redshift 2 = by applying a standard FoF 
analysis to the DM particle set, with linking length 0.15 times the 
mean inter-particle distance. After the FoF identification, the cen- 
tre of the halo was set to the position of its most bound particle. 
A spherical overdensity criterion was then applied to determine the 
virial radius, R200 of each halo. In this manner, 117 halos were 
identified within the simulated volume, of which 105 have virial 
mass A/200 greater than 10^^ /i^^Mq. To save computing time, 
the subsequent analysis was performed on a random subsample of 
roughly 2 million particles selected from the whole sample of 480'^ 
particles. We deliberately chose to analyse DM particles rather than 
galaxies in the simulations, as we do not trust the galaxy distribu- 
tion, with a homogeneous core found in these and other hydrody- 
namic cosmological si mulations, as they conflict with the observa- 
tions dSaro et al.|[20o3) . Moreover, the particle data provide much 
better statistics. 

A sample of 93 non-b imodal mock clusters were e x tracte d 
from these simulations by iMamon, Biviano, & Murant3 ( l2010h . 
who placed their observer at 90 Mpc of each mock cluster 
and added the peculiar motions to the distance with the Hubble 
constant value (Ho = 100 km s^^ Mpc^^) corresponding to the 
units of positions {h~^ k pc). This yields fa irly regular clusters that 
should resemble those of lYang et al. I( l2007l) (in the latter case, their 
high threshold causes multimodal clusters to be fragmented into 
several). We repeated the same extraction, this time with by plac- 
ing the observer at an infinite distance of each cluster (cylindrical 
projection). 



with projected radius and absolute LOS velocity. The generally in- 
creasing radial trend of /gores is clearly modulated by Iiilos |: at 
all projected radii, the value of /gores is higher (lower) for high 

(low) |llLOs|- 



4 MODELLING THE VELOCITY MODULATION 

In this section, we present simple models to quantify the observed 
fraction of GORES by combining theoretical predictions with the 
suite of galaxy properties analysed in ij3] 

We build two simple models of the frequency of GORES, one 
in which their frequency only depends on the class of galaxies, i.e. 
whether they are part of the virialised, infall or backsplash popu- 
lation, and second, where the frequency of GORES depends upon 
class and the distance from the centre of the cluster. We define the 
three classes of galaxies in radial phase space (i.e. [r, Vr], where r 
is the 3D radius and Vr the radial velocity) using the z = Q out- 
put of a cosmological simulation. We also use this simulation to 
project these quantities onto the observable projected phase space, 
[R, |wlos|], where R is the projected radius. 



4.1 Cosmological simulation 

We use the high-resolution cosm ological hydrodynamical simula- 
tion run bv lBorgani et alj ( [20041) . To summarise, the simulation as- 
sumes a cosmological model with fio = 0.3, = 0.7, Q,h = 0.039, 
/i = 0.7, and erg = 0.8. The box size is L = 192 h^^ Mpc. The simu- 
lation uses 480'^ DM particles and (initially) as many gas particles, 
for a DM particle mass of 4.62x10® h~^MQ. The softening length 
was set to 22.5 h'^ comoving kpc until redshift z — 2 and fixed 



4.2 Classification schemes 

Ideally, one should define the virial, inf all and backsplas h classes 
by following the orbits of the particles ( iGill et al.ll200i) . but this 
task is beyond the scope of the present article. Instead, we define 
the classes using the present-day radial phase space distribution of 
particles. Particles that are outside the virial radius, and travelling at 
velocities lower than the critical velocity, are falling into the clus- 
ter for the first time, while particles outside the virial radius with 
high er velocities must have crossed the cluster at least once. 

ISanchis, Lokas, & MamonI 12004) have analysed the distribu- 
tion of DM particles in the radial phase space (r, Vr) &X z = Q, for 
a cluster-mass halo in a dissipationless cosmological n-body simu- 
lation (2 56^ particles in a bo x of 150 Mpc width). Their fig. 1 
(see also lMamon et alj|2004b suggests that the separation of parti- 
cles belonging to the infall, backsplash and virialised populations 
can be written with a critical radial velocity delimiting the infall 
regime from the other two classes: 

Fig. |9] shows the distribution of DM particles in the stack of 
93 clusters in radial ph ase space {r,Vr). Overplotted are the cuts 
of ISanchis et alj | |2004|) designed to separate the virialised, infall, 
and backsplash classes of particles (Eq. [2j. The particles within 
the virial radius with large absolute radial velocities have uncertain 
classes. 

We devise several Sc hemes to identify th e three classes of par- 
ticles (we refer to fig. 6 in lBertschingeilll985l for a schematic view 
of the trajectories of particles in radial phase space). Table [T] sum- 
marises these schemes. The rapidly infalling particles with r < 
might represent the low-end tail of the radial velocity distribution 
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Table 1. Schemes of particle classes from positions in radial phase space 



Scheme Region 

r/-Rv 


A 

< 1 


B 

< 1 

(2) 

''r.crit ^ "r.crit 


C 

< 1 

r,crit 


D 

> 1 

(2) 

r,crit ' j^iii' 





virial 


virial 


virial 


infall 


1 


virial 


infall 


infall 


infall 


2 


backsplash 


infall 


infall 


infall 


3 


infall 


infall 


infall 


infall - 


4 


backsplash 


backsplash 


backsplash 


infall 


5 


backsplash 


backsplash 


infall 


backsplash 



Note: The letters A to D correspond to the regions in Fig.|9] The symbol o 
means 'in between' . 



of the virial class, which we denote by Scheme 0, or might in- 
deed be part of the infall class (Schemes 1 to 3), or may even be 
mainly populated by particles on their second infall, i.e. backsplash 
class (Scheme 4). The rapidly expanding particles with r < rv 
might again represent the high-end tail of the radial velocity distri- 
bution of the virial class (Schemes and 1), or the backsplash class 
(Schemes 2 and 4), or the infall class (if the star forming properties 
of galaxies in this region of radial phase space reflect those of the 
rapidly infalling class within the virial radius. Scheme 3). Scheme 
corresponds to maximum virialisation while Scheme 3 corresponds 
to maximum infall. 

Although Fig.|9]is based upon a very different simulation (hy- 
drodynamical instead of dissipationless and with three times better 
mass resolution) and on a stacked halo instead of a single one, it 
shows that the critical radial velocity of Eq. |2] is fully adequate 
to distinguish between infalling and backsplash particles in the 
stacked mock cluster obtained from the cosmological hydrodynam- 
ical simulation analysed here. 

One may argue that the red solid line separating the infall and 
backsplash classes is too high and that we may be missing an im- 
portant fraction of backsplash particles at high r and low Vr. We 
thus created a maximum backsplash Scheme 5, where we consider 
the critical velocity 



-^ = -2.2 + 0.83 



(3) 



This is shown as the magenta dashed line in Fig.|9l particles with 

(2) 

Vr < v^. are infalling, other particles are backsplash or virial 
according to Scheme 4 applied with Ur.crit- 

The cosmological simulation provides us with a unique way of 
measuring the frequency of the three classes of particles (galaxies) 
in projected phase space. We translate velocities from Vv ~ Vioo 
to the observable cluster velocity dispersion, a^, noting that the 
stacked cluster has a velocity dispersion (limited to the aperture 
r < Rv) of 0.65 Vv, as measu red in the hydrodynamical sim- 
ulation by iMamon et all bOlOh . and within 5% of what is ex- 
pected for isotropic NF W models of clusters (Appendix A of 
liviauduit & Mamotj|2007h . 

Fig. [To] shows the fractions of the three classes in projected 
phase space. The virialised particles are necessarily at radii r<Rv, 
the infalling particles prefer large projected distances and high ve- 
locities, while the backsplash particles are mostly found just around 
the virial radius and have low absolute LOS velocities. This prefer- 
ence of backsplash galaxies for low |i'los| is an immediate conse- 



backsplash 




Figure 9. Radial phase space distribution of dark matter particles of the 
stack of 93 regu lar mock clusters from the cosmological hydrodynamical 
simulation of ;Borgani et alj )2004l) . The units of radius and radial velocity 
are the virial radius rioo and the circular velocity at that radius, respec- 
tively. The critical velocity separating infall from backsplash population 
(long diagonal red line) is given in Eq. [2] The short diagonal line repre- 
sents the negative critical velocity threshold (see i|4j. For clarity, only 1 out 
of 550 particles of the original simulation is plotted. Region A is virialised 
in Schemes and 1, infall in Scheme 3, and otherwise backsplash. Region B 
is virialised in Scheme 0, infall in Schemes 1 to 3, otherwise backsplash. Re- 
gion C is virialised in Scheme 0, backsplash in Scheme 4, otherwise infall. 
Region D is backsplash in Scheme 5, otherwise infall (see Table[T). 



qu ence of their low er radial velocities (SjTJ and has also been noted 
bv lGill etalj ( l2005l) . 

In the top panel of Fig.[TT] we show the contours for the frac- 
tion of particles that are backsplash, while the bottom panel pro- 
vides the fraction of backsplash particles versus projected radius 
in wide bins of absolute LOS velocity. The two plots of Fig. [TT] 
indicate that the fraction of backsplash particles reaches a maxi- 
mum of 54% just outside the virial radius at very low Iulos] for 
Scheme 3 (75% for Scheme 5). At high |ulos|, the fraction of 
backsplash particles is symmetric relative to the virial radius. How- 
ever, at low |t'LOs|, the fraction is skewed towards higher values 
beyond the virial radius (see, e.g. curve representing all particles 
with |?;los| < 3 cr^ in the bottom panel of Fig.fTTll. Indeed, at pro- 
jected radii smaller than the virial radius, the backsplash particles 
must be in the foreground or background of the virial sphere, and 
geometric effects typically cause their 3D velocities to be aligned 
with the line of sight, making it difficult to obtain small absolute 
LOS velocities. These trends are confirmed in Table |2] which lists 
the fractions of the three classes of particles in cells of projected 
phase space. 



4.3 Models 

4.3.1 Global deprojection 

We first directly deproject the observed fractions of GORES, 
regardless of the LOS velocity. Expressing the global (sum- 
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Table 2. The fraction of virialised (v) infall (i) and backsplash (b) particles of the stacked mock cluster in cells of projected phase space with Scheme 3 



r/i?v 
^LOS 




0-0.5 






0.5-1 






1-1.5 






1.5-2 




V 


1 


b 


V 


1 


b 


V 


1 


b 


V 


1 


b 


0-1 


0.89 


0.08 


0.03 


0.49 


0.33 


0.18 


0.00 


0.65 


0.35 


0.00 


0.93 


0.07 


1-2 


0.83 


0.16 


0.02 


0.40 


0.48 


0.11 


0.00 


0.81 


0.19 


0.00 


0.97 


0.03 


2-3 


0.69 


0.31 


0.00 


0.16 


0.83 


0.01 


0.00 


0.99 


0.01 


0.00 


1.00 


0.00 
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Figure 10. Density in projected phase space of virial (green), infall (blue) 
and backsplash (black contours) particles of the stacked mock cluster with 
our preferred model (Scheme 3, see Fig. |9). The contours are logarithmi- 
cally spaced by a factor of 2.3, and the same set of contours is used for all 
three classes: the highest contours are on the lower left (virial), lower right 
(infall) and bottom (backsplash), where the latter two classes only reach the 
2nd highest contour of the virial class. 



ming over the velocity bins) observed fraction of GORES as 
5GOREs(i?) = NGOREs{R)/Ntot{R), we can deduce the surface 
densities, E(i?) = N(R)/{2n RdR) of GORES and of all galax- 
ies. Since the projected number profiles A''gores (R) and A^'tot (R) 
are fairly noisy, we fit a polynomial to log E vs. log R and extrapo- 
late the data beyond the last point by a power-law, whose slope is fit 
on the last three points. Fig. [T2] shows the surface density profiles 
and their polynomial fits. While the fits diverge dramatically at low 
R, they match remarkably well at high R (better than analogous 
fits of log E vs. R, not shown in Fig.[T2]for clarity). One can now 
perform Abel inversion to deduce the space densities 

u(r) = / (dE/di?)/\AR2372di? (4) 

Ti" Jr 

of GORES and all galaxies, and finally obtain the deprojected frac- 
tion of GORES by dividing the two space densities. Fig.|13|shows 
the resultant deprojected fraction of GORES, obtained with the or- 
ders 3 to 5 polynomial fits of log E vs. log R (Fig. I12l l. Fig. [T3] 
also shows a fit of Eq.|5]to the deprojected fraction of GORES 
obtained from the order 4 polynomial fit of log E vs. log R, using 



0.5 1 1.5 2 




0.5 1 1.5 2 



Figure 11. Upper panel: contours in projected phase space of the fraction 
of particles of the stacked mock cluster that contribute to the backsplash 
class (Scheme 3; see Fig.|9). Contours are linearly spaced from 0.05 to 0.5 
going downwards. Lower panel: fraction of backsplash particles in bins of 
absolute LOS velocity: 0-0.5 (top, black), 0.5-1.0 (red), 1.0-1.5 (green), 
1.5-2.0 (blue), and > 2.0 (t„ (bottom, magenta), and summed over all LOS 
velocities < 3 cr (thick brown). The dashed thick brown curx'e shows the 
maximum backsplash fraction (for IfLOsI < 3(t„) from Scheme 5 while 
the long-dashed thick grey line is the relation deduced from observations by 
iPimbblet. l .20lll) (see g3). 

/gores (r) = fo . (5) 

r + a 

We see that the simple, saturated linear model provided in Eq. [5] 
with asymptotic GORES fraction fo — 0.52 and quenching ra- 
dius for efficient star formation (hereafter quenching radius) a — 
1.26 Rv, is a decent representation of the deprojected profile for 
/gores- While this deprojection leads to an asymptotic value of 
/gores = 0.52, Eq.lHyields /gores = 0.44 at our typical mini- 
mum projected radius for field galaxies (6.7 Rv), very close to their 
observed fraction of 0.45. 

Note that although the constraints from potential SDSS fibre 
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5000 




Figure 12. Surface density of all galaxies (top) and of GORES (bottom). 
The symbols are the data (filled black circles for all galaxies and blue stars 
for GORES). The black, purple, and brown curves represent the polynomial 
fits of log E vs log R for orders 3, 4, and 5, respectively. The data were 
extrapolated beyond 2 with power-law fits to the outermost three data 
points. 



Table 3. Best-fit parameters for the 3 galaxy classes with the 6 schemes for 
the model of constant GORES fraction per class 



r < _Rv AND 



Scheme 



/v /i 



xi 






virial 


virial 


0.11 


0.33 


0.26 


4.6 


1 


infall 


virial 


0.10 


0.34 


0.24 


3.5 


2 


infall 


backsplash 


0.10 


0.33 


0.28 


2.9 


3 


infall 


infall 


0.10 


0.34 


0.23 


2.7 


4 


backsplash 


backsplash 


0.09 


0.32 


0.32 


3.0 


5 


backsplash* 


backsplash 


0.10 


0.33 


0.30 


2.9 



Notes: xi is the psr degree of freedom. The errors are typically 0.01 on 
/v and /i and 0.04 on /t,. In Scheme 5, the extreme negative velocities for 
r < Rv correspond to infall galaxies. 



collisions leads to underestimated surface density profiles at low 
projected radii, this bias in surface density should affect both the 
total and GORES surface densities in a similar (multiplicative) way 
(Eq.|4}. Hence the effect of fibre colfisions on the ratio of GORES 
to total space densities should be minor. 

We now refine the models accounting for the LOS velocity 
modulation of faoiiEsiR), using the dynamical classes of the 
galaxies (virial, infall, and backsplash). 




Figure 13. Deprojected fractions of GORES (using Eq. |4}. Thin cun'es 
show the results obtained by polynomial fits of the log surface density vs 
log R (orders 3, 4, and 5 in black, purple and brown, respectively). The 
thick green curve represents Eq.|5]with /o = 0.52 and a = 1.26 iJv, ob- 
tained by a fit to the order 4 polynomial fit of log S vs log R (purple 
curve), with r linearly spaced between 0.05 and 2 Rv 



4.3.2 Model 1: Constant recent starburst fraction per class 

Suppose that for each of the three classes of galaxies (virial, infall 
and backsplash), the fraction of galaxies which experienced a re- 
cent (< 1 - 3 Gyr) starburst (GORES) is independent of the 3D 
radius and radial velocity. Let these three fractions be /v, f\ and /b 
respectively (we drop the subscript 'GORES' from / for clarity). 

The fraction of GORES in a cell {R,Vz) of projected phase 
space can then be written as 

SCORES = '^faP{ce\Ri,Vj) , (6) 

a 

where, p{a\Ri,Vj) is the fraction of particles of class a within the 
cell {Ri,Vj) of projected phase space (as in Table 1, but with finer 
radial bins). 

One can then perform a fit between the observed and pre- 
dicted fractions of GORES gaoREsiRi , Vj ) over all cells [Ri ,Vj), 
using the conditional probabilities for each class, p{a\Ri, Vj), de- 
rived from the cosmological simulation. We assume binomial er- 
ror bars on the observed values of gGonEs{Ri,Vj) and that the 
p{a\Ri,Vj) derived from the cosmological simulation are perfectly 
known. Furthermore, we use the six Schemes for the three popu- 
lations, changing the class within the cluster virial sphere for the 
extreme velocities |tv| > |?'r,crit(f)| as explained in the caption of 
Fig.|9] Table [3] shows the best fits for each of these six population 
schemes. 

The best-fit predictions for the scheme with the lowest reduced 
(hereafter xl). Scheme 3 (where both the very positive and very 
negative radial velocity objects within the virial sphere are con- 
sidered to be part of the infall population), provide a reasonable 
match to the observed fractions of GORES in bins of projected 
phase space. But xt ~ 2.7 for this best-fit model, suggests that 
the model itself can be improved. 
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Figure 14. Same as Fig. [8] with best fit models overplotted: Top: tlie dashed 
and solid cui-ves sliow the best-fit Scheme 3 models (Tables[3]and|4j assum- 
ing no radial dependence (using Eq.|6) and the radially varying one of Eq|7] 
forcing Ov = Qi = f*b, using Eq. [s] respectively. Bottom: Best-fit radially 
varying model (Eq. [7) with free values of a^, a;, and a^, for Schemes 2 
(dashed) and 3 (solid). The schemes are defined in Tables[3]and|4] 



4.3.3 Model 2: Recent starburst fractions increasing with 
physical radius 

In our second model, we suppose that the fraction of GORES is 
no longer constant per class, but varies with physical radius r as in 
Eq.[5] varying the normalisation and scale for each class: 

r/iiv 



fa{r) = fa- 



il) 



r/ J?v + aa 

i.e. rising roughly linearly with radius for r < aa, and saturating 
to an asymptotic value fa at large radii. The predicted fractions of 
GORES is then 

5GOREs(-Ri,«j) = ^p(Qift,Uj)^/c,(rfc)(7(rfc|i?i,-!;j,a), 

(8) 

where q{k\Ri,Vj ,a) is the fraction of the particles of class a 
in the cell of projected phase space {Ri,Vj) that are in the fcth 
bin of physical radius. The physical radial bins are expressed 
as r-fc = Ri coshufe, where Uk is linearly spaced from to 
cosh" ^ r-niax / Ri , using r-max = 50 i?v . 

Table |4] shows the new set of fits. The first set of lines shows 
the fits when we force the same quenching radius for the three 
classes: Ov = ai = a^. Again, Scheme 3 provides the lowest x^, 
which is now 1.6 (in comparison with 2.7 for the / = constant 
model). So with the inclusion of just one extra parameter, the 



radially- varying model fits the data much better. This is confirmed 
with both the F-test and the AICc jAkaike|[T974l: iHurvich & Tsail 
Il989l) criterion (Table |4j. Note that cosmic variance on the finite 
number of clusters has been incorporated in the errors. The value of 
quenching radius a is low in Scheme 3, which indicates that it is not 
that far from the model with constant fraction of GORES per class. 
The top panel of Fig. [14] shows the observed fractions of GORES 
together with the best fits with Scheme 3 for both the constant frac- 
tion per class and radially- varying models with Ov — — a^^. 
Both models capture the initial rise of the GORES fraction with 
projected radius and subsequent saturation, as well as the velocity 
modulation, even though the fit is not excellent. 

In the next set of fits, shown in the bottom six rows of Ta- 
ble|4] we lift the assumption of av=ai — a^. The best-fit reduced 
values are now ^ 1.6 for Schemes 1-5, with an often signifi- 
cant improvement of the fit (despite the two extra parameters) over 
the radially-varying model with the same scale parameters (see the 
F test statistic in Table |4}. The i^-test indicates that the improve- 
ment is significant with the inclusion of two extra parameters (three 
scales instead of a unique one). Using the Akaike criterion, which 
reduces here to A AICc = AAIC — A^^, with associated prob- 
ability P — exp( — AAIC/2), we conclude that Scheme 3, which 
still leads to the lowest xt = 1-2, is a significantly better fit than 
Scheme (99.7% confidence) and marginally significantly better 
fit than Schemes 4 (93% confidence) and 5 (94% confidence). Al- 
though a decent fit to the data, this best-fit model still predicts 
too low values of /gores for high absolute LOS velocity galax- 
ies within half virial radius of the clusters. Interestingly, the best- 
fit model (Scheme 3) predicts a constant GORES fraction for the 
infall class galaxies (aj = 0), while the next best-fitting model 
(Scheme 2) predicts that both infall and backsplash populations 
have a constant GORES fraction (ai = Ob = 0). 



5 DISCUSSION 

This paper is a unique analytical effort to explore the variations in 
the spectral and physical properties of different dynamical classes 
of galaxies residing in similar environments projected on the ob- 
served sky (projected radius and absolute LOS velocity), and to 
interpret these differences. We draw a sample of ~ 20, 000 galax- 
ies (A4r ^ -20.5) found within 2R^ and IAulosI < Sa^ of 268 
galaxy clusters, and a similar sample of field galaxies for this pur- 
pose. All the galaxies are taken from the SDSS DR4 spectroscopic 
galaxy catalogue. 



5.1 Velocity segregation of stellar mass 

Three physical mechanisms might explain the excess of high mass 
galaxies in the cluster cores amongst the low velocity galaxies 
(Fig.|2](: cluster tides, two-body relaxation, and dynamical friction. 
We examine in turn each one of these. 

Ti dal effects from the c l uster gravitational field (predicted b 



iMerrittll 19831 : lMamorj[l993 : iGnedin. Hernquist. & Ostrikeil | 



clearly seen in cosmological simulations by Ghigna et alj 



eg by 
19991 



199^ 



and m easured through gravitational lensing by Nataraian et al.l 
l2009l) will be most effective on slowly moving galaxies, be- 
cause these have more time to feel the cluster tides. Even though 
stars are less affected by tides than the more extended dark mat- 
ter haloes, simulations of a live dwarf spiral galaxy orbiting 
a spiral galaxy ^100 times more mass ive, with fixed gravita- 
tional potential l lKlimentowski et al]l2009h or live A^-body system 
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Table 4. Best-fit parameters for the three galaxy classes with the four schemes for models of radially-varying GORES fraction 



r < i?v AND 



Scheme 


Vr < IV crit 


Vr > — fr.crit 


/v 


/i 


/b 




a; 




xl 


F 


Pf 


AAICc 


Paicc 









virial 


virial 


0.32 


0.38 


0.35 


0.58 


0.58 


0.58 


2.0 


46 





87 









1 


infall 


virial 


0.24 


0.37 


0.29 


0.36 


0.36 


0.36 


1.8 


31 





55 









2 


infall 


backsplash 


0.21 


0.36 


0.30 


0.28 


0.28 


0.28 


1.7 


24 





39 









3 


infall 


infall 


0.19 


0.36 


0.26 


0.22 


0.22 


0.22 


1.6 


24 





36 









4 


backsplash 


backsplash 


0.23 


0.35 


0.35 


0.35 


0.35 


0.35 


1.8 


24 





40 









5 


backsplash* 


backsplash 


0.22 


0.36 


0.34 


0.31 


0.31 


0.31 


1.8 


22 





37 





A1C-AIC(3) 


P 





virial 


virial 


0.41 


0.35 


1.00 


0.88 


0.27 


4.79 


2.0 


0.6 


0.56 


-3 




11.8 


0.003 


1 


infall 


virial 


0.34 


0.34 


0.36 


0.72 


0.00 


1.23 


1.6 


3.6 


0.04 


5.7 


0.06 


7.6 


0.02 


2 


infall 


backsplash 


0.30 


0.33 


0.21 


0.62 


0.00 


0.00 


1.5 


3.8 


0.03 


5.5 


0.06 


4.4 


0.11 


3 


infall 


infall 


0.29 


0.34 


0.39 


0.62 


0.00 


1.40 


1.2 


5.7 


0.01 


8.3 


0.02 






4 


backsplash 


backsplash 


0.26 


0.37 


0.28 


0.50 


0.55 


0.00 


1.6 


3.1 


0.06 


4.1 


0.13 


5.4 


0.07 


5 


backsplash* 


backsplash 


0.27 


0.35 


0.27 


0.53 


0.19 


0.00 


1.6 


2.8 


0.08 


3.4 


0.19 


5.7 


0.06 



Notes: The upper portion is for fits in which av = Qi = Qb is enforced, while the lower portion is for fits with fully free scale radii. The errors are typically 
0.01 on /v and /i and 0.04 on /b. In Scheme 5, the extreme negative velocities for r < Rv are infall. The best fitting scheme 3 is highlighted in bold. F and 
Pp represent the F test statistic that the fit with the scheme is significantly better than that of the previous model (/(r) with =0; =ab vs. / = cst for the 

upper panel and general /(r) vs. /(r) with av = = a^, for the lower panel) and the associated probabilities of a larger value o f F appearing by chance. 
AAlCc and PaiCc represent the same quantities with the Akaike criterion (AIC = x'^ + 2 k, for k Darameters. lAkaik3ll974h modified for finite size 
samples (AICc = + 2 k n/ (n — k — 1), for k parameters and n data points. iHurvich & Tsall989l) . The final two columns for the lower panel provide the 

Akaike statistic testing whether the fit for Scheme 3 is as good as Scheme 3 and its associated probability that a larger value of AICc occurs by chance. 



fcokas et alj|2010l) indicat e that, at each pericen tric passage, over 
35% I Klimentowski et aTl) to 67% fcokas et al.h of the stars of a 
spiral galaxy are tidally stripped. Tides will thus transform massive 
galaxies into low mass ones, and low mass galaxies into even lower 
mass ones. So tides cause a decrease in the fraction of high mass 
galaxies, which disagrees with the trend seen in Fig.|2] 

If two-body relaxation can lead to energy equipartition in the 
cluster cores, then the high mass galaxies will move slower (in 3D). 
So the fraction of high mass galaxies among the low 3D velocity 
ones should increase with decreasing cluster-centric radius. Hence, 
the fraction of high mass galaxies among the low absolute LOS 
velocity should also increase with decreasing radius. However, the 
same argument would lead us to predict that the fraction of high 
mass galaxies among the high velocity ones should decrease with 
decreasing radius, and instead Fig. |2] shows no such trend. There- 
fore, the trend in Fig. [2] cannot be due to two-body relaxation, 
unless some other process is increasing the fraction of high mass 
galaxies among the high velocity ones within the cluster core. 

The time required by dynamical friction to make the orbits of 
galaxies decay towards the cluster centre scales as /m, where v 
is the 3D velocity of the galaxy of mass m. Therefore the most mas- 
sive and slowly-moving galaxies should have their orbits decay the 
most, and be brought to even lower (projected) radius. So dynam- 
ical friction should boost the fraction of high-mass slow-moving 
galaxies in the cluster core. Moreover, the effects of dynamical fric- 
tion should be weak for the rapidly moving galaxies. This explains 
the lack of trend in the fraction of high mass galaxies with projected 
radius for the galaxies with large absolute LOS velocities. 

For dynamical friction (and tides), the high-velocity galaxies 
act little, if one assumes a population in virial equilibrium. Now, the 
infalling galaxies also tend to have high absolute radial velocities 
(Fig,|9]l, hence high absolute LOS velocities (e.g., Table[2}, in fact 
even higher than the virialised population. So the same argument 
applies to the infalling galaxies; their velocities are too high to be 
seriously affected by dynamical friction or tides. 

We note that the two-body relaxation timescale is independent 



of galaxy mass m, while the dynamical friction timescale varies 
as 1/m. It is easy to show that the ratio of two-body relaxation 
time to dynamical friction time scales as (m) /m. Thus, the lack of 
high-velocity low mass galaxies in the inner bin of projected radii 
in Fig.|2l in contrast with the excess expected from two-body relax- 
ation, is not inconsistent with the excess of high-mass galaxies at 
low radii observed in Figure |2] and expected from dynamical fric- 
tion. 

Is dynamical friction efficient enough to increase the frac- 
tion of high mass galaxies towards the centre among the slowly 
moving galaxies? We adopt the dynamical friction timescale that 
Ijiang et al] j2008l) carefully calibrated with a cosmological hydro- 
dynamical simulation: 



Ui ~ 1.4 



M{r)/r 



ln[l + M{r)/m] Vc{r) 



(9) 



where m is the galaxy mass (with its tidally stripped DM halo) and 
Vc{r) = \/ GM{r) jr is the circular velocity at radius r. Eq.|9]in- 
corporates the effects of elongated orbits and tidal stripping. One 
would like to have an idea of the physical radius of each galaxy 
given its projected radius and absolute LOS velocity. Fig.|15|shows 
the distribution of physical radii given the projected radii in our 
adopted bins of absolute LOS velocity. The physical radii of the low 
absolute LOS velocity galaxies are close to their projected radii, 
so that these low-velocity galaxies can be used to deproject radial 
trends. But high absolute LOS velocity galaxies appear to system- 
atically lie very far from the virial sphere (typically 10 i?v)- 

Using the relation r ~ 1.3R from Fig. [TS] we can esti- 
mate the dynamical friction time as a function of projected radius. 
Fig.[T6]indicates that a galaxy (along with its dark halo) must have 
a total mass of at least 10^^'^^ M0, for the orbital decay time, at 
R — 0. 25 Ry, to be small er than the age of the Universe. In com- 
parison, lYang et alj ( |2009|, their table 4) found that in the SDSS, 
central group/cluster galaxies of mass 3 x 10^'^M q resides in halos 
of mass 10^^'^ Mq. Also, according to table 4 of lYang et al.L cen- 
tral haloes of mass 1O^^ *^M0 have corresponding galaxy masses 
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Figure 15. Physical radius versus projected radius for particles of stacked 
mock cluster (medians with eiTor bars for quartiles) in bins of absolute LOS 
velocity: 0-1 (t„ {black triangles), 1-2 Cj, (red squares), and 2-3 (t„ {blue 
pentagons). The green curve shows r = R. 
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Figure 16. Timescale for orbital decay by dynamical friction versus pro- 
jected radius for galaxies of log total mass (in solar masses, including 
the tidally stripped dark matter halo) of (from top to bottom, labelled on 
right) 12 {red), 12.1 (green), 12.2 (blue), 12.3 (cyatt), 12.4 (magenta), and 
12.5 (orange), orbiting an NFW cluster of scale radius 0.19 i?v and mass 
My = Mioo = 1.9 X IO^'^Mq (corresponding to the median mass of 
Mi8o = lO^''-^ M0 of the SDSS clusters studied here) within the virial 
radius iJv = Rioo = 1.5 Mpc. The timescale is estimated using Eq.|2] 
assuming that r = 1.3 J? as inferred from our cosmological simulation (see 
Fig. list . The age of the Universe is shown as the horizontal line. 



of IO^^-^^Mq, i.e. 45% above our galaxy mass threshold. We there- 
fore conclude that dynamical friction might indeed be just suf- 
ficiently effective to explain the trends in Fig. [2] However, the 
DM halo of a galaxy extends further out than its stellar compo- 
nent, so halos are tidally stripped in a more efficient manner (e.g., 
iKlimentowski et al.l200a) . This means that the ratio of total to stel- 
lar mass of cluster galaxies is lower than for field galaxies. One 
would therefore expect the actual halo masses to be lower t han in - 
ferred for the central galaxies, from table 4 o f lYang et all ( |2009|) . 
and thus the orbital decay times by dynamical friction to be some- 
what longer. 

5.2 The star-forming galaxies 

The distributions of specific star formation rate, broad-band colour, 
and amplitude of the 4000 A break show statistical differences 
with absolute LOS velocity at almost all radii out to i? < 1.5 i?v 
(Figs. ID H and [6l Tables |A3l |A4] and |A6] respectively). How- 
ever, the field component shows a bimodal distribution in EW(H(5) 
(Fig. Is} and D,i4000 (Fig. [6ll around 2 A and 1.5 respectively. 
This is in agreement with the bimodality of these two diagnos- 
tics when plotted in terms of ste llar mass, luminosity, or con- 
centration jKauffmann et al.ll2003l) . We adopt these thresholds to 
define a galaxy as a recent starburst galaxy (GORES) if it has 
EW(H5) > 2 A and D„4000 < 1.5. 

As expected, the fraction of GORES drops within the virial ra- 
dius J?v for all velocity classes (Fig.[8j. But the fraction of GORES 
is modulated by the absolute LOS velocity: at a given projected ra- 
dius, galaxies with higher |?;los| have higher fractions of GORES. 
For the high absolute LOS velocity galaxies, the fraction of star- 
forming galaxies declines rapidly but only well within the core re- 
gion (< 0.5 i?v). This once again strengthens our interpretation that 
the high-velocity galaxies are the youngest members of the core, 
having recently fallen in. 

A curious feature occurs in the radial variation of both high- 
mass (Fig.|2ll and GORES fractions (Fig.[8j: the galaxies with high 
absolute LOS velocity show a marginally significant peak (dip) in 
the fraction of GORES (high mass galaxies) at \.5 R^. If this fea- 
ture is not a statistical fluke, it may be reflecting the discontinu- 
ity seen in SFR of cluster galaxies at the physical radius where 
infalling and virialised galaxies meet. However, this is expected 
to occur closer to the virial radius and our simple models involv- 
ing the virial, infall and backsplash populations do not lead to 
any discontinuity of /gores at 1.5 i?v (Fig.ll4t. An explanation 
may require a model going beyond spheri cal symmetry, e.g., with 
large- scale filaments feeding clusters (seelPorter & Ravchaudhurvl 
l2007l : |Porteretal.ll2008l : iMahaian et al.ll201l[) . For the benefit of 
the reader, we emphasise that the enhancement seen in the SFR of 
galaxies seems to occur in a very narrow region only. This obser- 
vation is very interesting and still being explored. It might prove 
useful to distinguish between the backsplash and infall galaxies. 

Given the M*-SSFR an ti-correlati on (e.g. 

Mamon. Parker. & Prous3l200ll : lNoeske et a l. 2007a; Damen et al.l 
20091) . the peak in the fraction of low mass galaxies might be 
directly linked to the enhancement of SFR in infalling galaxies on 
the outskirts of clusters before the environmental effects in clu sters 
predominate jPorter & Ravchaudhurvll2007l : |Porter et al.ll2008l) . 

While the projected radii of the low absolute LOS velocity 
particles in the stacked mock cluster underestimate only slightly 
their physical radii (Fig.llSK our cosmological simulation indicates 
that the high absolute LOS velocity particles with projected radii 
R > 0.8 Rv have very large physical radii (typically over 10 virial 
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radii). This implies that the observed fraction of GORES for the 
high absolute LOS velocity bin at i? > Rv should be the same as 
that of the field, but this is not what is observed (Fig.[8]l. One cannot 
resort to cosmic variance to explain this discrepancy, given the large 
number of clusters used in our analysis of the SDSS observations 
(268) and in the simulations (93), and we cannot find any other 
convincing explanation for this dichotomy. 



5.3 Identifying baclcsplash galaxies 

The analysis of the z — output of DM particles in our cosmo- 
logical hydrodynamical simulation allows us to distinguish back- 
splash particles, whose positions in radial phase space imply that 
they previously crossed their parent (mock) clusters from virialised 
particles (within the virial sphere) and the remaining infall particles 
(which are actually expanding away from their clusters when they 
are beyond the turnaround radius of ~ 3.5 7?v). Our analysis in pro- 
jected phase space indicates that the best place to detect backsplash 
particles is just beyond the virial radius and at very low velocities. 
With Scheme 3, the backsplash fraction reaches 54% (bottom panel 
of Fig.lllb for the lowest absolute LOS velocity bin. With the max- 
imum backsplash Scheme 5, it reaches 75% for the lowest |i)los1 
bin. 

IPimbbleJ ( l201ll) compared the distribution of SSFRs of galax- 
ies with 0.3 < |ulos|/o"u < 0.5, attributed to the infall class with 
the full distribution. This led him to predict that the fraction of 
backsplash galaxies decreases as /b = 0.61 — 0.052 R/Rv Such 
a trend (grey line in the bottom panel of Fig.ll 111 appears inconsis- 
tent with the backsplash fractions predicted from our cosmological 
simulation, even with our maximum backsplash Scheme 5. It is dif- 
ficult to understand how one would find that as many as half of all 
galaxi es at 2 viri al radii (the maximum projected radius in the anal- 
ysis of |Pimbble3) from clusters are backsplas h. Backsplash galaxi es 
typically only bounce out to 1-2 virial radii jMamon et al J 20041) . 

Yet, the high fraction of backsplash galaxies found by 
IPimbblej agrees w i th the 54±20% fraction of backsplash parti- 
cles that iGill et alj ( |2005|) found in the range of projected radii 
1 < R/Rv < 2 from their cosmological simulations. In our 
Schemes 0-4, we find that only 17%±1?^ of the particles in this 
range of projected radius (with |i;los| < 3 a„) are backsplash par- 
ticles. However, our backsplash fraction amongst DM particles is 
marginally consistent with the higher fraction found bv lGill et al.l 
Even in t he maximum backsplash Scheme 5 (partly motivated by 
IPimbble j * s high backsplash fraction, but only marginally consistent 
with the observed fraction of GORES), we find only 40±1% of the 
particles in the same range of projected radius are backs plash. This 
is apparently inconsistent with the fra ction deriv ed bv IPimbblell 
although fully consistent with that from Gill et alK simulations. 

Part of our disagreement with IPimbbleJ on the fraction of 
backsplash galaxies in the range 1 < R/R^ < 2 may be caused 
by the narrow range of IvlosI/o"!) used by him to select the infall 
galaxies. 



5.4 Tlie quencliing of star formation in 3D: inf ailing vs 
backsplash galaxies 

The SDSS spectroscopy provides the means to address the role of 
the dynamical classes (virial, infall and backsplash) in shaping the 

^ The errors in our fractions of backsplash particles are from 50 bootstraps 
on the 93 halos of the cosmological simulation. 



properties of galaxies by studying the spectral features as a function 
of both projected radius and absolute LOS velocity. 

Beyond the well-known radial dependence of galaxy proper- 
ties in clusters, our analysis highlights the velocity modulation of 
these trends for the parameters SFR/M*, EW(H(5) and the Dn4000. 
All trends with projected radius are amplified with increasing ab- 
solute LOS velocity. As can be seen in Fig.[T5] the effect of |ulos i 
effectively adding to the projected radius is not trivial. Neverthe- 
less, as qualitatively expected from the Hubble flow, in every bin of 
projected radius the higher absolute LOS velocities leads to higher 
physical radius. 

In ij3]we find observational evidence to support the statistical 
differences between the low velocity (|u|los ^ o"u) and the high ve- 
locity (2 < |wlos l/fi, ^5 3) galaxies on the outskirts of clusters (1- 
1.5 i?v). K-S tests between the different distributions in the clusters 
and in the field show statistically significant differences between 
the low and high velocity galaxies between 1-1.5-Rv in SFR/M*, 
EW(H(5) and the D„4000 (Tables lATlA6l l. 

Despite the noisiness of the variations with projected radius of 
the numbers of galaxies and of GORES, we have been able to de- 
project the fraction of GORES to show that it is well fit (Fig.ll 3t by 
a saturated linear model (Eq.[5]l. Such a model applied in turn to the 
three different classes of galaxy populations (virialised, infalling 
and backsplash) with the same quenching radius (qv = ai = a^), 
provides an adequate fit to the variation of GORES with projected 
radius in three bins of absolute LOS velocity (top four rows of Ta- 
ble |4] and upper panel of Fig. I14t. The model reproduces the rise 
with projected radius and the higher fractions for higher absolute 
LOS velocities. Similar models with free star formation quench- 
ing radii fit the data even better (bottom four rows of Tableland 
lower panel of Fig.ll4t. In particular, the radius-dependent model 
where galaxies within the virial sphere (r < _Rv) with radial veloc- 
ity \vr \ > liir.critl are cousidcrcd to be infalling, provides the best 
fit to the data (g31 Fig.|9]and Table|4j. 

In contrast, a simpler model, where the fraction of GORES 
is a constant that only depends on the class of galaxy population, 
does not provide an adequate fit (Table[3), even if it also reproduces 
these qualitative trends (Fig.ll4t. 

The ability of Scheme 3 to better fit the SDSS data than the 
other three Schemes suggests that the high absolute radial veloc- 
ity galaxies within the virial sphere, both infalling and outflowing, 
retain their fraction of GORES. However, our best-fit model with 
Scheme 2 is also an adequate fit, so the model with quenching at 
pericentre is also consistent with the data. 

In the constant GORES fraction model with Scheme 3, the 
GORES fraction in the backsplash class is half-way between the 
analogous values for the virial and infalling populations, suggest- 
ing that one passage through the cluster environment causes half 
the quenching of star formation. We now consider the deprojected 
models to assess how star formation is quenched at a unique radius 
where the three classes can be compared; the virial radius. 

Fig. Hz] shows how the best-fitting models predict the fraction 
of GORES as a function of the 3D radius. We first note that for the 
constant /gores per class model, the mean variation of /gores 
(top panel of Fig. I17t departs more strongly from the global trend 
found by the global deprojection ( i|4.3.U than do the mean varia- 
tions with the other models. This reinforces the view that the con- 
stant /gores per class model does not reproduce well our data. 

With the radially-dependant model of the GORES fraction 
with equal quenching scale radii, again with Scheme 3, the fraction 
of GORES in the backsplash population at r = i?v is intermediate 
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Figure 17. Fraction f(r) of GORES versus physical radius, for the virial 
{brown), infall (purple) and backsplash (blue-grey) populations, their mean 
(using the cosmological simulation to determine the radial variation of the 
densities of the three classes, black) and the global deprojection (green 
dash-dotted cun'e). Top: best fit constant (dashed lines) and radially in- 
creasing (Eq. |5) with Ov = di = CLh (solid curx'es), both with Scheme 3. 
Bottom: best fit radially increasing (Eq.|5) with free star formation quench- 
ing radii for Schemes 2 (dashed) and 3 (solid cuives). 



Table 5. Quenching at r = (models with free quenching radii) 



Scheme /q^resC^v) /S'oresC^v) /gores(^v) quenching 






0.22 


0.28 


0.17 


1.7 


1 


0.19 


0.34 


0.16 


1.3 


2 


0.18 


0.33 


0.21 


0.9 


3 


0.18 


0.34 


0.16 


1.1 


4 


0.17 


0.24 


0.28 


-0.7 


5 


0.18 


0.30 


0.27 


0.2 



Note: the last column lists Q = {fif^^Es - /gores)/(/gores " 
^GOHEs)' *hich measures the effectiveness of quenching of the back- 
splash population relative to the infall and virial ones. The best fitting 
scheme 3 is highlighted in bold. 



between those of the infalling and virialised populations, although 
closer to the virialised one (top panel of Fig. 117b. 

Table [5] displays the fraction of GORES at the surface of 
the virial sphere (r — R^) for the better-fitting radially- varying 
models with free quenching radii. The backsplash quenching fac- 
tor, Q = (/gores ~ /go'res)/(/gores ~ /gores)' should 
lie between (no quenching: backsplash and infall have the same 



GORES fraction at a given radius) and 1 (full quenching: back- 
splash and virial have same GORES fraction at a given r). The val- 
ues above unity imply that the fraction of GORES in the backsplash 
population lying on the virial sphere is lower than the correspond- 
ing fraction of the virial population at the same location. This leads 
to the unphysical result that on the virial sphere backsplash galaxies 
are more passive than the virialised ones. 

For our best-fitting Scheme 3, Q = 1.1 is just above unity, i.e. 
the backsplash and virialised populations on the virial sphere have 
very similar fractions of GOREqj, while for the second best fitting 
Scheme 2,Q = 0.9 (see bottom panel of Fig.fTTl). On the other hand, 
for Scheme 4, Q is much smaller than 0: the fraction of GORES 
for galaxies lying on the virial sphere is higher for the backsplash 
galaxies than for the infall population. This again appears to be an 
unphysical result, which leads us to disregard Scheme 4 (which fits 
the data worse than Scheme 3, but with only marginal significance; 
see Table|4}. This suggests that the backsplash galaxies are strongly 
quenched relative to the infalling ones, at least half the way from 
infalling to virialised galaxies at r = (for Scheme 5, which is 
only marginally consistent with the data) to nearly fully quenched 
to the level of virialised galaxies (Schemes 3 and 2). 

In contrast, usi ng hydrodynamical cosmological simulations 
of the Local Group, iKnebe et al.l j201lh find that the luminosity 
function of backsplash galaxies is similar to that of the infalling 
galaxies and that the total mass within the radius containing the 
outermost stars divided by the luminosity is higher for infalling 
galaxies than for backsplash galaxies. The first result suggests lit- 
tle quenching of star formation as dwarf galaxies cross through the 
virial sphere of the Milky Way, while the second result suggests 
inverse quenching, that is that backsplash galaxies have higher lu- 
minosity per unit mass. However, star formation in 10^^ Af© galax- 
ies falling into 10^^ Mq clusters should be similar to the quenching 
of star formation in dwarf galaxies of mass 10^° Mq falling into 
th e Milky Way virial sphere. What is more surprising is that fig. 5 
of lKnebeetatlfeOllb indicates that galaxies in the virialised class 
have higher luminosity for given internal velocity dispersion than 
galaxies in the infall class. Analysing the residuals of luminosity 
versus velocity dispersion relative to the mean global trend, we de- 
duce that the higher luminosity of virialised galaxies compared to 
infall galaxies of the same velocity dispersion is statistically sig- 
nificant (with 99.9% confidence using a K-S test, and by typically 
two magnitudes). The backsplash population on the other hand is 
intermediate but not inconsistent with either the virial or the in- 
fall populations. The higher luminosities of the virial population 
of simulated dwarfs appear to contradict the observed positive lu- 
minosity - radius correlation found for luminous virialised galaxies 
with clu ster-centric distance s above > 0. li?v from their parent rich 
clusters jAdami et"al]| 19981) . 

The quest for the physical mechanism responsible for this 
quenching is beyond the scope of this paper, be it tidal stripping, 
ram pressure stripping, harassment, starvation or pre-processing 
in groups. For instance, Ivan den Bosch et al] ( 1200 Sl) argue that the 
lack of dependence of galaxy colour on cluster mass is evidence 
against the dominance of ram pressure stripping and harassment. 

We note that in our two best-fitting models (Schemes 3 and 
2), the fraction of GORES in infalling galaxies is independent of 
their physical distance to the cluster. Our best-fit model requires a 
scheme (like our Scheme 3) where galaxies within the virial sphere 
with high outgoing velocities are equivalent to the infalling galax- 

® Forcing Q = 1 in our fit leads to xt only 0.003 higher for Scheme 3. 
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ies within the same volume (those with strong negative velocities). 
This together with the comparison of GORES fraction at the virial 
radius suggest that the star formation in galaxies is quenched on a 
timescale comparable to the time taken for a galaxy to bounce out 
of a cluster and reach its virial radius. This timescale corresponds 
to about 1 to 2 Gyr, which is very close to the lookback time to 
which our joint H(5-D„4000 GORES diagnostic is most sensitive 
( i]2.3t . Correcting for this lookback time, one concludes that star 
formation in a galaxy is nearly fully quenched in a single passage 
through the cluster. 

Of course, there will be a variety of orbital pericentres for the 
infalling galaxies, leading to different degrees of quenching. But 
the typical high elon gations found for orbits in ACDM simulations 
( iGhigna et al.lll998h suggests that orbits with large pericentres will 
be rare. 

We find that 17.9±0.4?{Q of the galaxies within the virial 
cylinder (projected radii R<Rv and \vlos | < 3 (t„) are GORES. 
This can be compared to the fraction of blue galaxies within the 
virial cylinder: we find that 20.3 ± 0.4% of the galaxies with 
R < Rv are 0.2 magnitude bluer than the Red Sequence (this frac- 
tion i^^irly_consistent with the ~78% of of SDSS cluster galaxies 
that I Yang et al.|[20o3 find to lie on the Red Sequence). So, within 
the virial cylinder, the fraction of GORES is roughly ScS±0.5% of 
the fraction of blue galaxies, although some blue galaxies are not 
GORES and some GORES are not blue jMahaian & Ravchaudhurvl 
l2009l) . 

In comparison, u sing the same cosm ological simulation as 
in the present article, iMamon et al. I J2O10h recently showed that 
23 ±1% of galaxies wit hin the virial c ylinder are outside the virial 
sphere. As discussed bv lMamon et al.L the match between fraction 
of blue galaxies and fraction of cluster interlopers is probably a 
coincidence since some galaxies within the cluster sphere must be 
blue and some interlopers must be Red Sequence galaxies in pro- 
jected groups. 

With our models, we can estimate the fraction of GORES 
within the virial sphere using the average of the mean fractions 
/(r) (Fig.|17t weighted by times the galaxy number density pro- 
file, v{r) measured in our cosmological simulation. Results given 
in Table |6] indicate that for our three best-fitting models (see Ta- 
bleQ, the fraction of GORES within the virial sphere is 13 ± 1%. 
We then predict the fraction of blue galaxies within the virial sphere 
to be 13/0.88 = 15 ±2%. 

Now, while most galaxies within the virial sphere are part of 
the virial population, some are infalling (Schemes 1 to 3) or back- 
splash (Scheme 2). Yet, we can proceed even further and estimate 
the fraction of GORES among the virialised class, again using our 
cosmological simulation to perform the calibration. As seen in Ta- 
ble |6l the fraction of GORES among the virialised class is not 
zero, but typically 11±1%. We also predict that the fraction of blue 
galaxies among the virialised population is 11/0.88 = 12.5 ± 2%. 
Moreover, Table [6] indicates that the fraction of GORES among 
the R < 2 Rv infall galaxies is 34 ± 1% while among the back- 
splash galaxies the fraction of GORES is typically 20 ± 4% (for 
Schemes 3, 2 and 1) or perhaps as high as 27 ± 4% (Schemes 4 and 
5). 

The presence of GORES among the virialised class might be 
caused by those low metallicity galaxies for which our H5-D,i4000 
estimator of recent efficient star formation has a lookback time as 



^ The errors here and below on the fractions of GORES are binomial and 
do not incorporate cosmic variance. 



Table 6. Fractions of Galaxies with Ongoing or Recent Efficient Star For- 
mation 



Model 


Scheme 


Virial 


Virial 


Virial 


Infall 


Backsplash 






cylinder 


sphere 


class 


class 


class 


Range 






r<_Rv 


(r<Rv) 


R<2Rv 


R<2Rv 


/=cst 


3 


0.17 


0.12 


0.10 


0.34 


0.23 




3 


0.17 


0.13 


0.11 


0.33 


0.22 


f{r) 


1 


0.17 


0.13 


0.12 


0.34 


0.19 


fir) 


2 


0.17 


0.13 


0.11 


0.33 


0.21 


fir) 


3 


0.17 


0.13 


0.11 


0.34 


0.19 


fir) 


4 


0.17 


0.13 


0.11 


0.32 


0.26 


fir) 


5 


0.17 


0.13 


0.11 


0.34 


0.27 


SDSS 




0.176 











Notes: The last six models are our best-fitting ones (xl < 1.6). The best 
fitting scheme 3 is highlighted in bold. Note that the fractions of GORES 
within the virial cylinder slightly underestimate the observed fraction, but 
ai'e consistent with it within the 1% errors. 

long as 3 Gyr ( i]2.3l l. Alternatively, these GORES might be the con- 
sequence of star formation triggered by rapid flybys and rare ongo- 
ing direct (satellite-satellite) mergers in clusters. 

Finally, our statistics lead us to conclude that a fraction 0.13 x 
(1 - 0.23)/0.18 = 0.56 of the GORES within the virial cylin- 
der are within the virial sphere, so that 44±2% of GORES within 
the virial cylinder are outside the virial sphere, for all good-fitting 
Schemes. 

Note that the asymptotic behaviour of the best-fitt models dis- 
agrees with the field value by at least 10%. One may think that this 
discrepancy may be due to a different set of environment-dependent 
phenomenon become efficient in modulating galaxy properties at 
R> 2 Rv from the centres of clusters, or to the oversimplified mod- 
els we used to deproject the fraction of GORES. However, accord- 
ing to the cosmological simulation we analysed here, the high ab- 
solute LOS velocity galaxies atR > Rv ought to lie at similar very 
far distances from the cluster centres (12 virial radii, see Fig. |15ll 
as the field galaxies. Therefore, the discontinuity in the observed 
fraction of GORES between these galaxies and the field is puzzling 
and remains an open question. 



5.5 Epilogue 

In summary, our work shows that the galaxy properties in and 
around clusters ar e not simple functions of stellar mass and loca l 
environment (e.g. iHaines et alJboOTi : Ivon der Linden et al. I|20l3), 
but also of absolute LOS velocity. With these velocities, it is then 
statistically feasible to segregate the infalling and backsplash galax- 
ies at any cluster-centric radius, including cluster outskirts. This 
velocity modulation will be better seen in the properties of the rel- 
atively low-mass galaxies, because they are the first ones to be 
influenced by any c hanges in their immediate environment (see 
iMahaian et al.ll201ll for such an effect seen in the nearby Coma 
cluster). It would be worthwhile to extend our analysis using other 
indicators of very recent star formation to 1) confirm our conclu- 
sion that the quenching of star formation occurs in a single passage 
through the cluster, and 2) assess the fraction of galaxies with very 
recent star formation among the virialised population. The survival 
of low-mass galaxies and of star formation in galaxies of all mass, 
as they pass through the cluster core, lead to important issues that 
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should also be addressed with high-resolution cosmological hydro- 
dynamical simulations. 



ACKNOWLEDGEMENTS 

SM is supported by grants from ORSAS, UK, and the Univer- 
sity of Birmingham. We thank the anonymous referee for useful 
comments and references. GAM thanks G. Murante for providing 
the cosmological simulation, and A. Biviano, G. Kauffmann, M. 
Sarzi, T. J. Ponman, and V. Wild for useful conversations. SM is 
grateful to B. M. Poggianti for a useful discussion. Funding for 
the Sloan Digital Sky Survey (SDSS) has been provided by the 
Alfred P. Sloan Foundation, the Participating Institutions, the Na- 
tional Aeronautics and Space Administration, the National Science 
Foundation, the U.S. Department of Energy, the Japanese Mon- 
bukagakusho, and the Max Planck Society. The SDSS Web site is 
|http://wvyw.sdss.org7| 



REFERENCES 

Adami C., Biviano A., Mazure A., 1998, A&A, 331, 439 
Adelman-McCarthy J. K., the SDSS collaboration , 2006, ApJS, 
162, 38 

Akaike H., 1974, IEEE Trans. Automatic Control, 19, 716 
Balogh M. et al., 2004a, MNRAS, 348, 1355 
Balogh M. L., Baldry I. K., Nichol R., Miller C, Bower R., Glaze- 
brook K., 2004b, ApJ, 615, L 101 
Balogh M. L., Navarro J. F, Morris S. L., 2000, ApJ, 540, 113 
Bekki K, 1999, ApJ, 510, L15 
Bertschinger E., 1985, ApJS, 58, 39 

Biviano A., Girardi M., Giuricin G., Mardirossian F, Mezzetti M., 
1992, ApJ, 396, 35 

Biviano A., Katgert P., Mazure A., Moles M., den Hartog R., 
Perea J., Focardi P, 1997, A&A, 321, 84 

Blanton M. R. et al., 2003, ApJ, 592, 819 

Blanton M. R. et al., 2005, AJ, 129, 2562 

Borgani S. et al., 2004, MNRAS, 348, 1078 

Bothun G. D., Schombert J. M., 1990, ApJ, 360, 436 

Brinchmann J., Chariot S., White S. D. M., Tremonti C, Kauff- 
mann G., Heckman T, Brinkmann J., 2004, MNRAS, 351, 1151 

Bruzual G., A., Chariot S., 1993, ApJ, 405, 538 

Chariot S., Fall S. M., 2000, ApJ, 539, 718 

Chariot S., Longhetti M., 2001, MNRAS, 323, 887 

Chincarini G., Rood H. J., 1977, ApJ, 214, 351 

Damen M., Labbe I., Franx M., van Dokkum P. G., Taylor E. N., 
GawiserE. J., 2009, ApJ, 690, 937 

de Vaucouleurs G., 1961, ApJS, 6, 213 

Dressier A., 1980, ApJ, 236, 351 

Ghigna S., Moore B., Governato F, Lake G., Quinn T, Stadel J., 

1998, MNRAS, 300, 146 
Gill S. P D., Knebe A., Gibson B. K., 2005, MNRAS, 356, 1327 
Girardi M., Rigoni E., Mardirossian R, Mezzetti M., 2003, A&A, 

406, 403 

Gnedin O. Y., Hemquist L., Ostriker J. P, 1999, ApJ, 514, 109 
Gunn J. E., Gott J. R., 1972, ApJ, 176, 1 

Haines C. P., Gargiulo A., La Barbera F, Mercurio A., Merluzzi 

P, Busarello G., 2007, MNRAS, 381, 7 
Haines C. P., La Barbera F, Mercurio A., Merluzzi P., Busarello 

G., 2006, ApJ, 647, L21 
Helou G., Salpeter E. E., Krumm N., 1979, ApJ, 228, LI 



Hurvich C. M., Tsai C.-L., 1989, Biometrika, 76, 297 

Jiang C. Y, Jing Y. P, Faltenbacher A., Lin W. R, Li C, 2008, 

ApJ, 675, 1095 
Kauffmann G. et al., 2003, MNRAS, 341, 33 
Klimentowski J., Lokas E. L., Kazantzidis S., Mayer L., Mamon 

G. A., 2009, MNRAS, 397, 2015 
Knebe A., Libeskind N. I., Knollmann S. R., Martinez- Vaquero 

L. A., Yepes G., Gottlober S., Hoffman Y, 201 1, MNRAS, 412, 

529 

Larson R. B., Tinsley B. M., Caldwell C. N., 1980, ApJ, 237, 692 
Lokas E. L., Kazantzidis S., Majewski S. R., Law D. R., Mayer 

L., Frinchaboy P M., 2010, ApJ, 725, 1516 
Mahajan S., Haines C. P, Raychaudhury S., 2010, MNRAS, 404, 

1745 

Mahajan S., Haines C. P, Raychaudhury S., 201 1, MNRAS, 412, 
1098 

Mahajan S., Raychaudhury S., 2009, MNRAS, 400, 687 
Mahajan S., Raychaudhury S., Pimbblet K. B., 201 1, MNRAS, to 

be submitted 
Mamon G. A., 1992, ApJ, 401, L3 

Mamon G. A., Biviano A., Murante G., 2010, A&A, 520, A30 
Mamon G. A., Parker Q. A., Proust D., 2001, Publications of the 

Astronomical Society of Australia, 18, 232 
Mamon G. A., Sanchis T, Salvador-Sole E., Solanes J. M., 2004, 

A&A, 414, 445 
Mauduit J.-C, Mamon G. A., 2007, A&A, 475, 169 
Merritt D., 1983, ApJ, 264, 24 

Mohr J. J., Geller M. J., Fabricant D. G., Wegner G., Thorstensen 

J., Richstone D. O., 1996, ApJ, 470, 724 
Moore B., Katz N., Lake G., Dressier A., Oemler A., 1996, Nat, 

379, 613 

Moss C, Dickens R. J., 1977, MNRAS, 178, 701 

Natarajan P., Kneib J., Small I., Treu T, Ellis R., Moran S., 

Limousin M., Czoske O., 2009, ApJ, 693, 970 
Navarro J. F, Frenk C. S., White S. D. M., 1996, ApJ, 462, 563 
Noeske K. G. et al., 2007a, ApJ, 660, L47 
Noeske K. G. et al., 2007b, ApJ, 660, L43 

Oemler A., Dressier A., Kelson D., Rigby J., Poggianti B. M., 

Fritz J., Morrison G., Small I., 2009, ApJ, 693, 152 
Pimbblet K. A., 201 1, MNRAS, 411, 2637 
Pimbblet K. A., Small I., Edge A. C, O'Hely E., Couch W. J., 

Zabludoff A. I., 2006, MNRAS, 366, 645 
Porter S. C, Raychaudhury S., 2007, MNRAS, 375, 1409 
Porter S. C, Raychaudhury S., Pimbblet K. A., Drinkwater M. J., 

2008, MNRAS, 388, 1152 
Rines K, Geller M. J., Kurtz M. J., Diaferio A., 2003, AJ, 126, 

2152 

Rines K, Geller M. J., Kurtz M. J., Diaferio A., 2005, AJ, 130, 
1482 

Sanchis T, Lokas E. L., Mamon G. A., 2004, MNRAS, 347, 1 198 

Sandage A., Tammann G. A., 1976, ApJ, 210, 7 

Saro A., Borgani S., Tornatore L., Dolag K., Murante G., Biviano 

A., Calura F, Chariot S., 2006, MNRAS, 373, 397 
Sodre L., Jr., Capelato H. V., Steiner J. E., Mazure A., 1989, AJ, 

97, 1279 
Struble M. F, 1979, AJ, 84, 27 

van den Bosch F. C, Aquino D., Yang X., Mo H. J., Pasquali A., 
Mcintosh D. H., Weinmann S. M., Kang X., 2008, MNRAS, 387, 
79 

Vollmer B., 2009, A&A, 502, 427 

von der Linden A., Wild V., Kauffmann G., White S. D. M., Wein- 
mann S., 2010, MNRAS, 404, 1231 



20 Mahajan, Mamon and Raychaudhury 



Worthey G., Ottaviani D. L., 1997, ApJS, 111, 377 
Yang X., Mo H. J., van den Bosch F. C, 2008, ApJ, 676, 248 
Yang X., Mo H. J., van den Bosch F. C, 2009, ApJ, 695, 900 
Yang X., Mo H. J., van den Bosch F. C, Pasquali A., Li C, Barden 
M., 2007, ApJ, 671, 153 



Velocity modulation of galaxy properties 21 



Table Al. K-S test probabilities for the null hypothesis that pairs of subsamples have distributions of A/* aiising from the same pai'ent distribution. Columns 
v\, v2, and v3 coiTespond to IfLOsl/fu = 0^1. 1^2, and 2-3, respectively. Statistically significant differences between subsample pairs are highlighted in 
bold. 





R/Rv 5*0.5 


Q.5<R/R^ <:i.O 


1.0<R/i?v ^1.5 


1.5<i?/i?v ^2.0 


field 




vl v2 v3 


vl 


v2 


v3 


vl 


v2 


v3 


vl 


v2 


v3 






I 2 3 


4 


5 


6 


1 


8 


9 


10 


11 


12 




1 


l.OE+00 7.3E-08 2.SE-09 


4.2E-20 


2.9E-20 


1.6E-09 


1.2E-12 


1.8E-09 


4.1E-08 


2.7E-11 


1.4E-10 


8.5E-06 


O.OE+00 


2 


l.OE+00 1.7E-02 


3.1E-03 


6.2E-06 


9.9E-04 


5.9E-03 


6.1E-02 


1.2E-03 


2.1E-02 


5.5E-02 


6.5E-02 


2.6E-13 


3 


l.OE+00 


1.7E-01 


3.6E-01 


3.6E-01 


4.3E-02 


1.2E-01 


2.6E-01 


3.6E-02 


3.2E-01 


7.3E-01 


2.6E-01 


4 




l.OE+00 


2.8E-02 


1.3E-01 


2.1E-01 


7.6E-01 


2.3E-01 


1.4E-01 


8.8E-01 


6.7E-01 


5.6E-07 


5 






l.OE+00 


6.5E-01 


2.5E-02 


8.8E-02 


8.9E-01 


3.3E-03 


3.2E-02 


6.4E-01 


1.3E-01 


6 








l.OE+00 


3.8E-02 


9.9E-02 


7.8E-01 


2.0E-02 


7.9E-02 


2.6E-01 


1.8E-01 


7 










l.OE+00 


7.9E-01 


l.OE-01 


9.7E-01 


5.8E-01 


6.7E-01 


2.6E-04 


8 












l.OE+00 


1.5E-01 


8.4E-01 


7.9E-01 


8.3E-01 


5.4E-03 


9 














l.OE+00 


4.3E-02 


9.9E-02 


5.2E-01 


4.3E-01 


10 
















l.OE+00 


4.9E-01 


7.0E-01 


5.4E-06 


11 


















l.OE+00 


9.1E-01 


3.7E-03 


12 




















l.OE+00 


3.8E-01 


Table A2. Same as Table|AT] but with BCGs excluded. 




R/R^ ^0.5 


0.5<R/R^ ^1.0 


1.0<R/R^ ^1.5 


1.5<R/R^ s:2.0 


field 




vl v2 v3 


vl 


v2 


v3 


vl 


v2 


v3 


vl 


v2 


v3 






1 2 3 


4 


5 


6 


1 


8 


9 


10 


11 


12 




1 


l.OE+00 5.2E-08 6.6E-10 


l.SE-20 


1.2E-21 


1.2E-09 


6.9E-15 


9.0E-10 


6.SE-09 


2.0E-12 


9.5E-11 


3.9E-06 


O.OE+00 


2 


l.OE+00 9.2E-03 


1.5E-03 


1.7E-06 


2.3E-04 


8.9E-04 


3.6E-02 


5.1E-04 


4.0E-03 


3.8E-02 


3.8E-02 


3.9E-13 


3 


l.OE+00 


2.0E-01 


4.3E-01 
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9.8E-02 
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4.4E-02 


2.2E-01 


7.4E-01 


3.4E-01 
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l.OE+00 
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6.3E-02 
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8.5E-01 


6.5E-01 


2.6E-06 
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4.6E-02 
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5.4E-03 
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6.7E-01 


5.8E-01 


6 
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6.7E-02 


7.7E-02 


7.3E-01 


3.4E-02 


4.1E-02 
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1.5E-01 
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8.9E-01 


8.1E-01 


7.8E-01 


4.5E-03 
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1.3E-01 


8.9E-01 


7.3E-01 


8.7E-01 


1.7E-02 
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5.2E-02 


6.9E-02 


4.5E-01 


1.9E-01 
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4.9E-01 


6.8E-01 


5.8E-05 
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l.OE+00 
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8.2E-03 
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l.OE+00 


5.3E-01 



APPENDIX A: KOLMOGOROV-SMIRNOV STATISTICS 

In this appendix, we present tables of probabilities, from the Kolmogorov-Smirnov test, that various sub-samples of different physical 
quantities arise from the same parent distributions (see §|3j. 
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Table A3. Same as Table lAll but for the difference between the (g — r)''-^ colour and the best-fitted red sequence (RS). 



H/i?v ^0.5 ().5<R/Rv <1.0 1.0<_R/_Rv ^1.5 1.5<_R/_Rv ^2.0 field 





vl v2 


v3 


i;l 


v2 


v3 


v\ 


v2 


vi 


v\ 


v2 


v3 






1 2 


3 


4 


5 


6 


1 


8 


9 


10 


11 


12 




1 


I.OE+OO 1.2E-01 


l.OE-04 


5.9E-30 


3.4E-29 


1.3E-10 


O.OE+00 


O.OE+00 


2.0E-29 


O.OE+00 


O.OE+00 


9.0E-27 


O.OE+00 


2 


l.OE+00 


1.7E-02 


1.5E-16 


6.1E-17 


1.8E-07 


l.OE-41 


7.3E-41 


6.2E-23 


1.4E-45 


8.4E-36 


7.6E-20 


O.OE+00 


3 




l.OE+00 


5.8E-05 


4.6E-06 


1.6E-02 


1.2E-13 


4.SE-16 


1.8E-12 


3.9E-15 


2.5E-13 


2.4E-08 


O.OE+00 


4 






l.OE+00 


5.6E-02 


3.3E-01 


1.3E-11 


1.2E-14 


3.8E-10 


2.2E-15 


2.2E-11 


l.lE-07 


O.OE+00 


5 








l.OE+00 


8.1E-01 


l.lE-04 


1.6E-07 


1.3E-06 


5.4E-06 


1.8E-05 


l.OE-04 


O.OE+00 


6 










l.OE+00 


2.4E-03 


4.1E-05 


6.1E-05 


8.0E-04 


9.9E-04 


5.5E-03 


4.9E-29 


7 












l.OE+00 


l.lE-01 


1.8E-02 


6.1E-01 


2.7E-01 


l.OE-01 


O.OE+00 


8 














l.OE+00 


3.3E-01 


1.3E-01 


6.5E-01 


3.6E-01 


1.7E-23 


9 
















l.OE+00 


5.0E-02 


1.5E-01 


1.5E-01 


3.0E-06 


10 


















I.OE+OO 


8.1E-01 


4.8E-01 


O.OE+00 


11 




















l.OE+00 


5.6E-01 


1.4E-31 


12 






















l.OE+00 


1.9E-16 



Table A4. Same as Table|AT] but for SFR/A/* . 



R/i?v s;0.5 0.5<i?/flv <1.0 1.0<iJ/i?,v ^1.5 1.5<ii:/i?v ^2.0 field 





vl v2 


v3 


i;l 


v2 


v3 


vl 


v2 


vl, 


vl 


v2 


vi 






1 2 


3 


4 


5 


6 


1 


8 


9 


10 


11 


12 




1 


l.OE+00 5.8E-04 


l.OE-04 


l.OE-26 


8.1E-28 


4.7E-12 


O.OE+00 


5.6E-45 


l.lE-22 


O.OE+00 


1.8E-40 


1.5E-21 


O.OE+00 


2 


l.OE+00 


1.7E-02 


l.OE-08 


6.2E-13 


1.3E-06 


4.1E-25 


l.lE-25 


1.2E-16 


2.1E-33 


1.7E-21 


1.2E-14 


O.OE+00 


3 




l.OE+00 


3.4E-03 


1.4E-02 


3.1E-02 


1.5E-06 


l.lE-07 


7.SE-08 


l.OE-07 


2.5E-06 


3.4E-06 


1.6E-21 


4 






l.OE+00 


6.1E-02 


8.4E-02 


4.1E-07 


1.6E-09 


4.SE-09 


2.3E-11 


3.2E-08 


7.9E-07 


O.OE+00 


5 








l.OE+00 


4.7E-01 


6.2E-04 


l.lE-04 


6.1E-06 


3.0E-05 


3.4E-04 


1.8E-04 


1.4E-26 


6 










l.OE+00 


9.1E-03 


2.3E-02 


6.3E-03 


4.1E-03 


1.6E-02 


l.lE-02 


5.5E-07 


7 












l.OE+00 


1.5E-01 


2.5E-03 


3.0E-01 


7.3E-02 


2.9E-02 


1.9E-17 


8 














l.OE+00 


9.7E-02 


4.6E-01 


5.2E-01 


4.5E-01 


2.9E-0S 


9 
















l.OE+00 


1.7E-02 


8.9E-02 


7.5E-01 


l.lE-01 


10 


















l.OE+00 


4.3E-01 


1.9E-01 


2.7E-14 


11 




















l.OE+00 


6.0E-01 


1.8E-08 
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Table AS. Same as TableEH but for H(5. 



R/i?v <0.5 0.5<R/R^ ^1.0 1.0<R/R^ s£1.5 1.5<R/R^ s£2.0 field 





til v2 


v3 


t>l 


v2 


v3 


t>l 


v2 


v3 


1)1 


v2 


v3 






I 2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




1 


l.OE+00 1.7E-06 


7.0E-09 


O.OE+00 


O.OE+00 


6.1E-22 


O.OE+00 


O.OE+00 


8.9E-37 


O.OE+00 


O.OE+00 


l.lE-27 


O.OE+00 


2 


l.OE+00 


6.6E-04 


1.4E-18 


2.6E-24 


9.6E-14 


3.7E-38 


7.7E-39 


4.4E-25 


1.5E-39 


6.2E-3S 


1.6E-17 


O.OE+00 


3 




l.OE+00 


6.0E-03 


8.5E-05 


2.8E-03 


7.3E-07 


2.8E-10 


2.0E-09 


2.7E-06 


9.4E-08 


1.2E-04 


4.4E-36 


4 






l.OE+00 


2.4E-03 


5.8E-03 


7.2E-07 


1.9E-11 


4.4E-11 


1.3E-08 


3.6E-09 


1.4E-04 


O.OE+00 


5 








l.OE+00 


7.8E-01 


2.0E-01 


2.4E-04 


2.7E-06 


5.5E-02 


1.3E-02 


7.8E-02 


1.6E-42 


6 










l.OE+00 


2.1E-01 


1.6E-02 


6.4E-04 


2.2E-01 


7.7E-02 


1.7E-01 


l.OE-12 


7 












l.OE+00 


3.6E-02 


3.6E-04 


l.OE-01 


3.1E-01 


5.7E-01 


6.1E-42 


8 














l.OE+00 


1.6E-01 


1.8E-02 


3.2E-01 


4.3E-01 


5.4E-13 


9 
















l.OE+00 


9.9E-04 


l.OE-02 


6.6E-02 


1.3E-02 


10 


















l.OE+00 


4.1E-01 


8.6E-01 


l.OE-43 


11 




















l.OE+00 


9.8E-01 


9.1E-19 


12 






















l.OE+00 


7.3E-09 



Velocity modulation of galaxy properties 

A6. Same as Table|AT] but for D„4000. 



R/R^ ^0.5 0.5<R/R-^ ^1.0 1.0<_R/_Rv =£1-5 1.5<R/R-^ sS2.0 field 





111 v2 




vl 


v2 


v3 


t'l 


v2 


v3 


vl 


v2 


v3 






1 2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 




1 


l.OE+00 3.5E-09 


5.1E-09 


O.OE+00 


O.OE+00 


9.8E-15 


O.OE+00 


O.OE+00 


8.8E-35 


O.OE+00 


O.OE+00 


1.2E-27 


O.OE+00 


2 


l.OE+00 


l.lE-04 


3.1E-23 


1.9E-23 


8.4E-09 


4.7E-40 


6.2E-38 


1.9E-24 


1.4E-45 


1.3E-35 


5.2E-18 


O.OE+00 


3 




l.OE+00 


8.1E-03 


7.2E-05 


4.7E-02 


7.4E-07 


9.2E-09 


5.7E-09 


9.0E-07 


3.7E-08 


8.2E-05 


1.2E-34 


4 






l.OE+00 


2.2E-02 


l.OE-01 


1.2E-06 


4.8E-09 


2.5E-07 


1.6E-11 


5.6E-09 


1.9E-05 


O.OE+00 


5 








l.OE+00 


2.8E-01 


6.3E-02 


8.2E-04 


9.8E-05 


l.OE-03 


5.8E-03 


1.7E-02 


O.OE+00 


6 










l.OE+00 


1.4E-02 


5.6E-04 


5.1E-05 


1.7E-02 


5.2E-03 


9.5E-02 


l.lE-17 


7 












l.OE+00 


7.1E-02 


7.5E-03 


l.lE-01 


6.9E-02 


1.8E-01 


O.OE+00 


8 














l.OE+00 


2.1E-01 


1.2E-01 


9.4E-01 


l.OE-01 


7.4E-17 


9 
















l.OE+00 


2.1E-02 


l.lE-01 


4.9E-02 


1.3E-03 


10 


















l.OE+00 


3.5E-01 


5.9E-01 


2.5E-38 


11 




















l.OE+00 


1.2E-01 


1.5E-21 


12 






















l.OE+00 


4.1E-10 



